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aerospace  industry  an  excellent  desig^analysis  tool  for  aircraft  such  as  the  Hi^  Speed  Civil 
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control  technology. 


TECHNICAL  ABSTRACT  (SUMMARY  REPORT) 

The  objective  of  the  Phase  I  STTR  contract  is  to  sufficiently  develop  the  Aeroservoelastic 
(ASE)  discipline  and  to  the  grounds  for  its  eventoal  integration  into  ASTROS  (Automated 
Structural  Optimization  System).  The  ZONA  Team  attempts  to  achieve  two  major  teclmical 
goals  in  order  to  accomplish  this  objective. 

First,  a  Unified  S-Domain  Aerodynamics  (USDA)  package,  as  derived  from  the  UAIC  (the 
Unified  Aerodynamic  Influence  Coefficient)  approach  of  four  major  ZONA  aerodynamic 
software  (called  the  ZAERO  Module)  covering  the  complete  flight  Mach  number  regime,  have 
been  developed.  A  feasibility  study  has  been  conduct^  for  th^  wing  planforms,  namely,  i) 
the  NASP  Demonstrator  Model,  ii)  the  AGARD  Standard  445.6  Wing  (solid  and  weakened 
planforms)  and  iii)  the  Model^  F-16  Wing.  Emphasis  was  placed  on  the  feasibility  of 
employing  the  developed  USDA,  via  Karpel's  Minimum  State  Techitique  (MIST),  for  transonic 
and  hypersonic  flutter  investigation  of  these  wing  planforms.  The  feasibility  study  was 
concluded  successfully  assuring  the  USDA's  capability  for  ASE  applications. 

Second,  a  Software  Design  Blueprint  for  the  ASE  Module  architecture  and  plans  for  its 
integration  into  ASTROS  are  provided.  Thus,  with  the  approved  Software  Design  Blueprint  for 
ASE  as  a  basis,  the  subsequent  integration  of  ASTROS/ASE  can  be  implemented.  Technically 
or  commercially,  the  resulting  integrated  ASTROS  code  will  acquire  a  competitive  edge  over 
all  existing  MDO  software. 

In  parallel  to  the  above  effort,  additional  tasks  were  conducted  to  further  substantiate  the 
present  teclmical  goals;  these  are: 

i)  Transoitic  unsteady  pressure  results  of  the  Lessing  Wing  and  the  LANN  Wing:  For  further 
validation  of  the  ZTAIC  code  (see  Appendix  A). 

ii)  Development  planning  of  an  Aerodynamic  Geometry  Module  (AGM)  that  will  provide 
ASTROS  with  a  universal  set  of  geometric  definitions  for  the  ZA^O  module. 

The  findings  of  the  Phase  I  studies  can  be  summarized  as  follows: 

On  The  ZAERO  Module: 

i)  The  unified  AIC  formulation  is  ideal  for  a  MDO  environment  such  as  the  one  provided  by 
ASTROS. 

ii)  The  unified  feature  of  ZAERO  for  all  Mach  numbers  will  surpass  existing  aerodynamic 
modules  including  those  presently  within  ASTROS. 

iii)  The  success  of  the  USDA  results  ensure  the  feasibility  of  ASTROS/ASE  applications. 

On  The  ASE  Blueprint: 

The  approach  and  the  overall  capability  of  the  ASE  in  preparation  for  future  integration  with 
ASTROS  has  been  defined  (see  Section  4.0).  Four  specific  scenarios  pertinent  to  the 
ASTROS/ASE  applications  clearly  elucidate  its  future  capability. 
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SECTION  1 


INTRODUCTION 

Current  performance  of  advanced  fighters  and  bombers  requires  high  maneuverability,  agility 
and  stealth  under  a  wide  range  of  flight  conditions.  Their  design  goal  tends  to  achieve  at  more 
flexible  and  less  inherently  stable  aircraft  throughout  the  complete 
(subsonic/transonic/supersonic)  flight  regime. 

One  such  advanced  conceptual  design  is  the  Active  Aeroelastic  Wing  (AAW,  Refs  1,  2), 
currently  being  pursued  by  the  Air  Force  and  Rockwell/North  American.  A  recent  design 
concept  of  the  high  Speed  Civil  Transport  (HSCT)  also  places  particular  emphasis  on  the  fly-by¬ 
wire  control  of  the  akcraft  flexibility  and  its  relaxed  stability.  Clearly,  a  system  that  could 
handle  such  interactions  between  structures,  unsteady  aerodynamics  and  active  control  is 
warranted.  Such  a  system,  known  as  Aeroservoelasticity  (ASE),  is  a  Multidisciplinary 
Technology  and  has  been  in  rapid  progress  over  the  last  two  decades  (Refs  3-9). 

Active  Aeroelastic  Wing  (AAW)  technology  is  a  multi-disciplinary  synergistic  technology  in 
which  the  proper  handling  of  the  interaction  between  aerodynamics,  control  and  structures  would 
allow  a  thin,  high  aspect  ratio  wing  planform  to  be  aeroelastically-deformed  in  achieving 
optimum  performance.  Software  support  for  the  AAW  design,  such  as  ASTROS  with  an  ASE 
capability,  would  therefore  be  most  desirable. 

There  exist  two  major  developments  in  aeroelastic  codes  that  can  handle  ASE  to  a  certain  extent, 
namely,  the  ISAC  code  developed  by  NASA/Langley  (Ref  3)  and  MSC/NASTRAN  Version  68 
(Ref  10).  However,  neither  one  can  handle  the  design/analysis  ASE  problem  adequately  in  that 
the  ISAC  code  lacks  the  capability  in  optimization  and  design  and  MSC/NASTRAN's  ASE 
capability  is  not  yet  fully  developed.  In  the  last  decade,  ASTROS  has  made  rapid  progress  (Refs 
1 1-19)  and  has  demonstrated  its  outstanding  capability  in  many  areas  of  MDO  technology  except 
ASE.  Thus,  the  proposed  Phase  I  work  in  the  development  and  integration  of  the  ASE  into 
ASTROS,  once  accomplished,  will  become  a  cutting-edge  software  technology. 

In  Phase  I  of  S'lT'R  (AF95T(X)9),  ZONA  Technology,  Inc.  and  its  team  members,  Oklahoma 
University,  Universal  Analytics,  Inc.  and  Dr.  Mordechay  Karpel  (the  ZONA  Team)  have  been 
committed  to  the  ASE  development  in  ASTROS  and  have  worked  out  the  ASE  Software  Design 
Blueprint  for  the  ASE  module  architecture  and  its  integration  into  ASTROS.  In  parallel  to  this 
effort,  a  prototypical  ZONA  Aerodynamic  (ZAERO)  module  (Refs  19, 20),  which  would  readily 
generate  the  Unified  Aerodynamics  (UMC)  for  all  Mach  numbers  and  Unified  S-Domain 
Aerodynamics  (USDA)  for  ASE  applications,  has  been  developed  and  properly  interfaced  with 
ASTROS.  A  feasibility  study,  which  included  the  testing  of  USDA  as  a  tool  for  the  flutter  on 
three  wing  planforms  extending  to  the  Transonic  and  Hypersonic  flow  regimes,  was  conducted 
successfully.  From  these  studies,  it  was  demonstrated  that  the  ZAERO  module  is  by  far  a  more 
general  tool  than  the  current  Aero  Module  in  ASTROS.  Also,  it  has  become  apparent  that  the 
ZAERO-generated  USDA  is  an  excellent  tool  for  ASE  application,  valid  throughout  the 
complete  flight  regimes  (Ref  20). 
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SECTION  2 


THE  ZAERO  MODULE 
2.1  Overview  of  the  ZAERO  Module 

The  ZAERO  module  consists  of  four  major  unsteady  aerodynamic  codes  that  jointly  cover  the 
complete  domain  of  all  Mach  number  ranges,  namely  ZONA51U,  ZONA6,  ZONA7  and  ZTAIC. 
As  can  be  seen  in  Figme  2.1.1,  the  AERO  modules  currently  integrated  within  MSC/NASTRAN 
and  ASTROS  only  possess  purely  subsonic  and  supersonic  capabUities. 


MSC/ 

NASTRAN 

AERO 


NONE 

1 

ZDNA51 

1 

0 

5 

MACH  NUMBER 


ASTROS 

AERO 

MODULE 


NONE 


Figure  2.1.1  ZONA  Aerodynamic  Module. 


By  contrast,  the  ZAERO  module  serves  as  a  unified  aerodynamic  tool  which  provides  computed 
data  from  unsteady  pressures  to  GAFs  throughout  all  Mach  numbers  by  means  of  the  unified 
AIC  approach.  In  fact,  it  is  the  UAIC  of  the  ZAERO  Module  that  has  efficiently  provided  the  k- 
domain  solution ,  and,  hence,  the  s-domain  solution  for  subsequent  ASE  application. 

The  development  of  codes  in  the  ZAERO  module  have  been  the  major  endeavor  of  ZONA 
Technology  in  the  last  decade.  The  following  is  a  brief  account  of  the  capabilities  of  these  four 
computer  codes;  namely 

•  ZONA51U 

•  ZONA7IZONA6 

•  ZTAIC 
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•  Z0NA51 U:  Generates  Unified  Unsteady  Hypersonic! Supersonic  Aerodynamics  for  Lifting 
Surface  Systems  (Refs  21, 22). 


A  Unified  Supersonic/Hypersonic  Lifting  Surface  Method  has  been  developed  by  ZONA 
recently  (Refs  21,  22).  This  method  combines  the  Supersonic  Lifting  Siuface  Theory  (such  as 
ZONA51,  Ref  23)  with  a  nonlinear  correction  matrix  based  on  Donov  &  Linnell's  uniformly 
valid-higher-order  Hypersonic/Supersonic  scheme.  This  correction  matrix  takes  the  flow 
nonlinearity  as  well  as  the  flow  rotationality  due  to  oscillatory  shock  waves  into  account,  which 
covers  both  the  Mach  wave  and  Newtonian  limits.  For  aeroelastic  applications,  ZONA51U  has 
been  applied  to  various  wing  planforms  with  thickness  distributions  (e.g.  Rectangular  Wings 
with  Diamond/Wedge  profiles  and  a  70'  Delta  Wing,  see  Figs  2.1.2  and  2.1.3,  Refs  21,  22). 
Computed  results  are  found  to  agree  well  with  those  computed  by  Euler  methods;  flutter  results 
are  validated  with  measured  data.  It  is  found  that  ZONA51U  improves  substantially  over  the 
linear  theory  results  in  terms  of  pressures  and  stability  derivatives,  and  provides  more 
conservative  flutter  boundaries  due  to  the  thickness  effect.  Furthermore,  the  input  format  of 
ZONA51U  is  nearly  the  same  as  that  of  ZONA51  with  only  an  additional  input  card  on  the  Wing 
Profile  Slope.  The  CPU  time  for  ZONA51U  is  also  comparable  to  that  of  ZONA51. 
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M 


(C) 

Figure  2.1.2  Damp-in-pitch  Derivatives  of  a  Rectangular  Wing  with  Diamond  Airfoil  Section  (h/c=0.5)  (a)  Cm^ 
vs.  Semi-Wedge  Angle  at  M=2.0  (b)  Cm^  vs.  Semi-Wedge  Angle  at  M=10.0  and  (c)  Ca/j  vs.  Mach 
Number  at  Semi-Wedge  Angle  of  15*. 


Figure  2.1.3  Flutter  Boundary  of  a  70*  Delta  Wing  with  6%  Biconvex  Airfoil  Section. 
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•  ZONA7IZONA6:  Generates  Unsteady  Supersonic! Subsonic  Aerodynamics  for  Aircraft 

Configurations  with  External  Stores  (Refs  26, 27, 28, 29). 


Prior  to  1990,  all  unsteady  aerodynamic  methods  for  aeroelastic  computations  were  based  on 
lifting-smface  models  (e.g.  DLM,  Ref  30).  The  aerodynamic  effects  due  to  the  presence  of 
bodies  and  due  to  the  wing-body  interference  are  largely  ignored.  However,  coupled  external- 
store  wing  flutter,  a  problem  that  is  of  frequent  concern  to  the  Air  Force,  cannot  be  resolved  by 
the  lifting  surface  modeling  alone.  To  demonstrate  this  effect.  Fig  2.1.4  shows  the  stability 
derivatives  of  a  NACA  wing-body  combination  based  on  wing-only  and  wing-body  analysis 
using  ZONA7.  Fig  2.1.5  presents  the  spanwise  imsteady  forces  and  moments  of  a  MLR  wing  (F- 
5  wing)  with  underwing  fin-missile  and  pylon.  It  is  seen  that,  in  both  cases,  the  discrepancies 
between  the  wing-only  results  and  the  wing-body  results  are  substantial. 

ZONA6  is  the  subsonic  counterpart  of  ZONA7  except  that  it  includes  the  important  body- wake 
effect  for  fuselage  and  stores.  It  should  also  be  noted  that  ZONA6’s  lifting  surface  option 
(referred  to  as  ZONA61)  has  the  same  order  in  paneling  scheme  as  that  of  ZONA51  and 
therefore  is  more  robust  scheme  than  that  of  DLM  (Ref  22). 

Fig  2.1.6  presents  the  out-of-phase  pressures  on  two  spanwise  stations  on  a  70°  Delta  Wing.  It  is 
seen  that,  using  typical  panel  cuts,  DLM  breaks  down  at  M  =  0.8  and  k  =  0.5. 

Fig  2.1.7  presents  the  unsteady  pressure  along  the  underwing  store  of  a  NLR  Wing-Tiptank- 
Pylon-Store  configuration.  It  is  clearly  seen  that  large  discrepancies  exists  between  the  body- 
only  and  the  wing-body  results. 
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PRESENT  MCTHOO  EXPERIMENT: 


A- 2  TRIANGULAR  WING  (NACA  0003-63  SECTIONS). 

(a) 


PRESENT  METHOD  EXPERIMENT^ 
— WING  +  BODY  OR«1.I8*IO^ 

-  WING  ONLY  o  R*  1.89  x  10 

-  BODY  ONLY 

-I.Er 


0.4' - ' - ' - ' - ' - ’ 

1.0  1.2  1.4  1.6  1.8  2.0 

MACH  NUMBER.  M 
(b) 


♦  WING  ♦BODY  Q  R»  1.18x10 

- V/ING  ONLY  O  R»  1.89x10 

- BODY  ONLY  O  R= 3.77x10 

- FORCE  TEST  DAT;. 

NACA  RM  A50K24: 


(C) 


Figure  2.1.4  Stability  Derivatives  of  a  NACA  Wing-Body  Combination  AR=2.0  Computed  by  Z0NA7  (a)  Sketch 
of  a  Wing-Body  Combination:  Delta  Wing  with  Centered  Body  of  Revolution  0))  Damping-in-Pitch 
Moment  Coefficients  with  Pitching  Axis  at  0.35c  (c)  Moment-Curve  Slopes  with  Pitching  Axis  at 
0.35c. 


PRESENT  EXP  (UNSTEADY)  EXP  (QUASI- STEADY  ) 

CLEAN  WING  -  - •  •  • - 

WING  ♦  PYLON  ♦  STORE  -  - +  +  _  — ■ - • ■ - 


Figure  2.1.5  Unsteady  Spanwise  Normal  Force  and  Pitching  Moment  for  the  Clean  F-5  Wing  and  the  Underwing 
Store  Configuration  at  M=1.35  and  Reduced  Frequency  k=0.1. 
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10x10 


Station  @ 


40x10 


Figure  2.1 .6  Out-of-Phase  Pressures  on  Two  Spanwise  Stations:  70'  Delta  Wing  (M=0.8,  k=0.5,  x^sO.Sc). 


(a)  (b) 

Figure  2.1.7  Unsteady  Pressure  Distributions  Along  the  Store  of  NLR  Wing-Tiptank-Pylon  Store  Configuration 
at  (a)  6=157.5  deg  and  (b)  6=292.5  deg  M=0.45,  a=0  deg.,  k=0.305  and  x^=0.15Cr. 
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•  ZTAIC:  Generates  Unsteady  Transonic  Aerodynamics  for  Lifting  Surface  Systems 
(Refs  32, 33). 


Since  1985,  ZONA  has  been  following  up  on  the  development  of  the  Transonic  Strip  (TES) 
Method  originally  supported  by  NAVAIR/ONR  (Refs  31,  32)  for  unsteady  flow  computations  of 
arbitrary  wing  planforms.  The  TES  method  consists  of  two  consecutive  steps  added  to  a  given 
nonlinear  Transonic  Small  Disturbance  Code  such  as  ZTRAN  (Ref  34),  namely  the  chordwise 
mean  flow  correction  and  the  spanwise  phase  correction.  Based  on  the  TES  concept,  ZONA’s 
Transonic  Aerodynamic  Influence  Coefficient  (ZTAIC)  method  is  developed  to  fully  automate 
the  computation  procedure  resulting  in  a  AIC  matrix  (Ref  33).  The  computation  procedure 
requires  direct  pressure  input  from  other  computed  or  measured  data.  Otherwise,  it  does  not 
require  airfoil  shape  or  grid  generation  for  a  given  planform.  Meanwhile,  all  the  mean-flow 
shock  jumps  are  properly  included  in  the  resulting  unsteady  aerodynamics  through  the  AIC 
formulation.  Results  computed  by  ZTAIC  have  been  validated  with  existing  results  for  a  number 
of  wing  planforms  with  various  aspect  ratios.  These  include  (shown  here):  the  Lessing  Wing  at 
M=0.9  (Ref  35,  Fig  2.1.8);  the  LANN  Wing  at  M=0.82  (Ref  36,  Fig  2.1.9);  and  the  Norhtrop  F-5 
Wing  with  control  surface  (Ref  37,  Fig  2.1.10).  More  validated  cases  as  a  part  of  the  cases 
studied  in  Phase  I  will  be  shown  in  following  sections.  ZONA’s  essential  improvement  lies  in 
the  ZTAIC  code  development,  which  not  only  can  provide  compatible  transonic  AICs  in  accord 
with  ASTROS'  format  but  also  achieves  considerable  savings  in  computer  time  and  reasonable 
accuracy.  Furthermore,  ZTAIC  has  a  user  oriented  input  format  which  is  fully  compatible  with 
that  of  DLM. 
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NASA  TND-344  /  H.C.  Lessing  &  J.L.  Troutman  /  M=0.24-1.3 
AR  ~  3.0, 5%  Parabolic  Arc 


Seeben  0  @  0 

n  -  O.S  0.7  0.9 


Section  1 


(a) 

Section  2 


Section  3 


1  - ZTAIC 

1  ®  Exp{  1  at  Run)  j 

1  ■  Exp{2nd  Run)  / 
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^  •  •  V  ' 

^ -  -  _ _  o 
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Section  1 
11  =  0.5 


-  ZTAIC 
Exptlsi  Run) 
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0.2  0.4  o.«  0.8 

Jt/c 


Section  2 
ti  =  0.7 


I  0  0.2  0.4  O.S  0.8 

x/c 


Section  3 
ti  =  0.9 


>  0  0.2 


0.4  0.8  O.S 

X/C 


Figure  2.1.8  The  Lessing  Wing  (a)  Configuration  (b)  Magnitude  of  the  First  Bending  (c)  Phase  Angle  (in  degrees) 
of  the  First  Bending,  M=0.9,  k=0.13. 
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2.2  Feasibility  Study  of  ZAERO/IJSDA  Methodology 


In  Phase  I,  the  proven  feasibility  study  of  the  USDA  focused  mainly  on  the  unified 
hypersonic/supersonic  and  transonic  flow  regimes  for  proper  verification  with  the  existing 
measured  data  and  CFD  results.  Three  major  wing  planforms  were  considered: 


•  the  NASP  Demonstrator  Model  (Ref  38) 

•  the  AGARD  Standard  445.6  Wing  (Solid  and  Weakened  Planforms)  (Ref  39) 

•  the  Modeled  F-16  Wing  (Ref  40) 


In  all  of  the  figures,  the  ZONA51U  and  ZTAIC  codes  were  used  to  generate  the  k-domain 
aerodynamics,  first  by  means  of  ZONA's  UAIC  formulation.  The  USDA  for  each  case  is  then 
obtained  through  Karpel's  expedient  s-domain  fit,  called  MIST  (the  Minimum  State  Technique, 
see  Eq  (4.6.2),  Refs  41, 42),  on  the  GAF’s  of  the  k-domain  aerodynamics.  Finally,  the  s-domain 
flutter  results  are  presented  in  Root-Locus  plots  and  were  compared  with  those  generated  by  the 
V-g  method. 
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•  NASP  Demonstrator  Wing 


Fig  2.2.1  shows  the  physical  dimensions  and  the  aerodynamic  modeling  of  this  planform.  Notice 
that  the  wing  cross-section  is  a  hexagon  profile. 

Fig  2.2.2  presents  the  results  of  the  Generalized  Aerod5mamic  Forces  (GAF's)  for  modes  3  and  4 
in  the  s-domain  at  three  hypersonic  Mach  numbers,  namely  M^o  =  5.0, 10.0  and  15.0. 

In  obtaining  these  results,  MIST  requires  only  five  aerodynamic  approximated  roots  covering  a 
range  of  reduced  frequency  from  k  =  0  to  k=  0.75.  It  is  seen  that  the  correlation  between  the  k- 
domain  results  and  that  of  the  USD  A  are  excellent. 

Fig  2.2.3  presents  the  Root-Locus  plot  at  these  three  Mach  numbers  for  8  modes.  For  cases  at 
Moo  =  5.0, 10.0  and  15.0,  mode  4  is  found  to  flutter  at  dynamic  pressures  at  130, 184  and  213 
psi,  respectively. 

Extensive  comparison  of  various  computed  results  has  been  carried  out  in  the  present  NASP 
wing  case.  In  Table  2.2.1,  the  flutter  results  of  ZONA51U  using  the  V-g  method  (k-domain)  and 
using  USDA  (s-domain)  are  compared  with  those  of  Piston  theory  and  the  NASA  CFD  codes 
(QSCFD2d  and  3d,  Ref  43).  Several  observations  can  be  made  with  respect  to  the  collected  data: 

-  Good  agreement  is  found  between  the  flutter  q’s  of  the  V-g  method  and  the  s-domain  method 
based  on  ZONA5  lU.  Flutter  frequencies  are  also  found  to  agree  well. 

-  Piston  theory  results  follow  the  general  trend  predicted  by  ZONA5  lU,  except  that  the  former 
results  tend  to  become  overestimated  at  high  Mach  numbers. 

-  Substantial  discrepancies  are  found  between  the  NASA/CFD  results  and  the  ZONA51U 
results:  ZONA51U  predicts  more  conservative  flutter  spe^s  at  higher  altimdes,  whereas  the 
NASA/CFD  codes  predict  less  conservative  speeds  at  lower  altitudes.  In  fact,  ZONA51U 
gives  the  most  conservative  answers  for  all  cases  in  this  smdy.  Our  confidence  on  these 
flutter  results  stems  from  our  thorough  verification/validation  of  the  ZONA51U  code  (Refs 
21, 22).  However,  further  validation  with  test  data  is  warranted  for  this  particular  smdy. 
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IIVIAGCS)  (rad/s) 


(a)  M=5.0  (b)  M=10.0  (c)  M=15.0 

Figure  2.2.3  Root-Locus  Plots  of  NASP  Demonstrator  Model  at  M=5.0, 10.0  and  15.0. 

(a)  8  Modes,  Qf  =  130  psi,  (Of  =  82  r/s  on  Mode  4 

(b)  8  Modes,  Qf  =  184  psi,  (Of = 7 1  r/s  on  Mode  4 

(c)  8  Modes,  Qf  =  213  psi,  (Of = 68  r/s  on  Mode  4 


Aerodynamic 

5 

Mach  Number 

10 

15 

Method 

qf 

psi 

r/s 

h 

(Kft) 

% 

psi 

C0£ 

r/s 

h 

(Kft) 

qf 

psi 

©f 

r/s 

h 

(Kft) 

Piston  Theory 

129 

78 

18 

184 

78 

42 

250 

72 

51 

QSCFD  2d 

169 

80 

11 

331 

81 

28 

982 

224 

22 

QSCFD  3d 

— 

— 

! 

330 

82 

28 

981 

224 

22 

ZONA51U 
-  V-g  Method 

129 

82 

18 

175 

73 

42 

206 

70 

55 

-  S'Domain 

130 

82 

18 

184 

71 

41 

213 

68 

54 

h  =  Approximate  matchpoint  altitude. 


Table  2.2.1  Comparison  of  Hutter  Dynamic  Pressures  and  Frequencies  of  NASP  Demonstrator  Model  at  M=5.0, 

10.0  and  15.0. 
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•  The  AGARD  Standard  445.6  Wing 


The  AGARD  Standard  445.6  wing  planform  has  two  structural  models:  the  solid  wing  and  the 
weakened  wing.  This  is  an  ideal  case  for  testing  ZTAICs  AIC  capability,  since  the  aerodynamic 
shapes  of  the  two  structural  models  are  identical  and,  therefore,  the  AIC  matrix  remains  the 
same. 

Fig  2.2.4  presents  the  flutter  results  of  the  weakened  wing.  At  a  subsonic  Mach  number,  say  M^o 
=  0.678,  tihe  ZTAIC  result  is  in  perfect  agreement  with  that  of  ZONA6,  as  expected.  At  other 
supercritical  transonic  Mach  numbers,  say  Mo©  =  0.95,  ZTAIC  predicts  a  pronounced  transonic 
dip  which  is  comparable  to  results  predicted  by  the  CAP-TSD  code  (Ref  44). 

Fig  2.2.5  presents  the  flutter  results  for  the  solid  wing.  Since  the  mean  planform  remained 
unaltered  for  the  aerodynamics,  the  AICs  used  here  would  be  the  same  as  those  used  for  its 
weakened  counterpart.  Thus,  only  a  re-start  of  the  code  (using  the  previously  generated  AIC 
matrix)  was  required  for  all  remaining  cases.  This  amounted  to  only  one  minute  of  CPU  time  per 
case,  a  great  savings.  Again,  the  ZTAIC  results  for  these  cases  agree  well  with  the  CAP-TSD 
result  and  measured  data.  This  demonstrates  that  ZTAIC  is  ideally  suited  for  aeroelastic 
optimization,  particularly  in  the  ASTROS  MDO  environment.  Fig  2.2.6  presents  the  MIST  fit  of 
the  GAFs  for  the  445.6  weakened  wing  at  M©©  =  0.95.  Unlike  tiie  previous  cases,  some  slight 
discrepancy  is  found  in  the  GAF  (1,2),  which  might  be  caused  the  stringent  feature  of  the 
supercritici  mean  flow  at  this  Mach  number. 

Fig  2.2.7  presents  the  root-locus  plot  of  the  same  case  using  the  MIST  fit  results  of  the  above. 
The  flutter  speed  predicted  here  (Vy  =  1015  fps)  is  in  reasonable  agreement  with  the  one  predicted 
by  the  V-g  method  Wf  =  944  fps,  see  Table  of  Fig  2.2.4).  The  discrepancy  is  caused  by  slight 
discrepancies  in  the  MiST-fit  in  the  transonic  range  (e.g.  see  Q12  of  Fig  2.2.6). 
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Figure  2.2.4  Comparison  of  Flutter  Speed  and  Frequency  of  445.6  Weakened  Wing  at  M=0.678, 0.90  and  0.95. 
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Figure  2.2.5  Comparison  of  Flutter  Speed  and  Frequency  of  445.6  Solid  Wing  at  M=0.90  and  0.95. 
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Figure  2.2.6  Minimum-State  Approximation  of  Generalized  Aerodynamic  Forces  of  445.6  Weakened  Wing  at 
M=0.95  Computed  by  ZTAIC. 
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Figure  2.2.7  Root-Locus  Plots  of  445.6  Weakened  Wing  at  M=0.95  Using  S-Domain  Aerodynamics  Computed 
by  ZTAIC  and  ZONA6. 
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Modeled  F-16  Wing 


Fig  2.2.8  presents  the  structural  model  and  the  flutter  boundary  of  a  modeled  F-16  wing.  It  is 
seen  that  the  flutter  results  of  ZTAIC  are  in  good  agreement  with  those  of  XTRAN3  and 
CAPTSD,  particularly  at  Moo  =  0.95.  It  should  be  noted  that  the  flutter  mechanism  changes  from 
Moo  =  0.9  to  Moo  =  0.95,  as  indicated  by  the  jump  in  the  flutter  frequency  (from  6.8  Hz  to  19 
Hz).  ZTAIC  concurs  with  the  other  two  codes  in  the  prediction  of  this  mechanism. 

Fig  2.2.9  presents  the  MIST-fit  of  the  s-domain  results  for  four  GAF's.  For  this  case,  MIST 
achieves  a  better  fit  than  for  the  case  of  the  445.6  wing.  Notice  that  the  dotted  lines  in  the 
background  is  the  fitted  subsonic  GAF's  using  ZONA6.  Substantial  differences  can  be  seen 
between  the  results  of  ZTAIC  and  ZONA6,  showing  strong  transonic  characteristics  for  this  wing 
at  Moo  =  0.95. 
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Flutter  Speed  and  Frequency  vs.  Mach  Number 
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-ZTAIC 


Figure  2.2.8  Flutter  Speeds  and  Frequencies  of  Modeled  F-16  Wing  Computed  by  ZTAIC,  CAPTSD  and 
XIPANS. 
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Figure  2.2.9  Minimum-State  Approximations  of  Generalized  Aerodynamic  Forces  of  Modeled  F-16  Wing  at 
M=0.95  Computed  by  ZTAIC  and  ZONA6. 
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SECTION  3 


AERODYNAMIC  GEOMETRY  MODULE  (AGM) 

3.1  The  Aerodynamic  Geometry  Module  of  ZAERO 

The  development  of  the  Aerodynamic  Geometry  Module  (AGM)  is  motivated  by  the  lack  of  a 
unified  input  format  for  various  aerodynamic  methods  which  would  compute  the  steady  and/or 
unsteady  aerodynamics  in  their  valid  Mach  number  ranges.  Specifically  different  flight 
conditions  require  different  computational  methods  and,  consequently,  require  different 
geometric  parameters  of  the  aerodynamic  configuration.  The  proposed  AGM  will  provide 
ASTROS'  AERO  module  with  a  universal  set  of  geometric  definitions  for  the  ZONA6,  ZONA7, 
ZONA51U  and  ZTAIC  codes  in  ZAERO  Module.  With  the  understanding  that  the  US  S AERO 
code,  currently  imbedded  in  ASTROS’  AERO  Module,  will  be  replaced  by  other  higher-order 
panel  methods,  the  AGM  will  also  have  the  capability  of  interfacing  with  other  higher-order 
panel  methods  such  as  QUADPAN,  PANAIR  and  ZONAIR  (Refs  24, 25). 

With  respect  to  the  geometric  definitions,  there  are  generally  three  types  of  panel  methods  (for 
both  steady  and  unsteady  codes): 

(i)  Lifting  Surface  Methods:  The  singularity  is  placed  on  the  mean-surface  of  the  lifting 
surface  such  as  the  Doublet  Lattice  Method  (DLM),  Constant  Pressure  Method  (CPM), 
ZTAIC,  ZONA51U,  ,and  the  Carmichael  flat  plate  method. 

(ii)  High  -Order  Surface  Panel  Methods:  The  singularity  is  placed  on  the  exact  surface  of 
the  configuration.  There  is  no  discrimination  between  body-like  components  and  wing¬ 
like  components.  These  methods  include  QUADPAN,  PANAIR  and  ZONAIR. 

(iii)  Hybrid  Methods:  The  singularity  is  placed  on  the  mean  surface  of  the  wing-like 
components  and  on  the  exact  surface  of  the  body-like  components.  These  methods 
include  ZONA6  and  ZONA7. 

Based  on  the  aforementioned,  two  types  of  input  will  be  incorporated  into  AGM:  the  "surface- 
panel"  input  for  higher-order  panel  methods  and  the  "thin-wing"  input  for  the  lifting  surface 
methods.  For  hybrid  methods,  no  additional  type  of  input  is  needed,  since  the  connectivity 
between  the  "thin-wing"  and  the  "surface-panel"  can  be  suitably  defined. 

Surface  Panel  Input 

The  most  convenient  way  to  define  surface  panels  for  modeling  an  arbitrary  configuration  is  the 
NASTRAN  type  of  input.  Five  new  bulk  data  entries  will  be  employed:  AGRID,  AQUAD, 
ATRIA,  AQDMEM  and  AB  AR. 

AGRID  defines  a  location  of  an  aerodynamic  grid.  AQUAD  and  ATRIA  define  a  quadrilateral 
and  triangular  surface  panels,  respectively.  AQDMEM  is  used  to  model  an  arbitrary  wake  panel. 
ABAR  defines  a  flat  wake  sheet  extending  from  the  wing  trailing  edge  or  the  end  of  a  truncated 
body  to  infinity.  Fig  (3.1.1)  demonstrates  the  use  of  these  new  bulk  data  entries  in  modeling  an 
arbitrary  wing-body-tail  configuration. 
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ABAR 


Figure  3.1.1  New  Bulk  Data  Entries  for  Aerod5Tiamic  Modeling  of  an  Arbitrary  Wing-Body-Tail  ’ 
Configuration. 

The  boundaty  condition  of  a  solid  body  in  inviscid  flow  is  satisfied  by  imposing  the  zero-normal- 
flow  condition  (no  flow  can  penetrate  into  the  body)  on  the  body  surface.  However,  this  is  not 
true  for  an  inlet  face  where  the  flow  is  allowed  to  penetrate  into  the  body.  To  define  this  kind  of 
boundary  condition,  all  AQUAD  and  ATRIA  cards  are  referred  to  as  a  set  of  "property  cards" 
defined  by  the  PAERO  bulk  data  en^.  The  PAERO  card  allows  the  user  to  specify  the 
percentage  of  the  mass  flow  ratio  that  is  allowed  to  penetrate  into  the  body.  The  value  of  this 
percentage  can  be  found  based  on  the  engine  operating  condition  (e.g.,  for  zero  mass  flow  ratio, 
the  zero-normal-flow  condition  is  imposed). 

The  user  may  often  be  interested  in  computing  aerodynamic  loads  on  a  group  of  components 
(e.g.  the  loads  on  an  external  store).  A  new  bulk  data  entry  "ACORD"  is  defined  which  specifies 
the  identification  number  of  a  local  coordinate  system  and  the  identification  number  of  a  group 
of  components. 

In  what  follows,  detailed  definitions  for  each  of  the  above  new  bulk  data  entries  are  given. 

•  Input  Data  Entry:  ACRID 


Description:  Defines  the  location  of  a  geometric  grid  point  of  the  aerodynamic  model. 


Format  and  Exaitple 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ACRID 

ID 

CORD 

X 

Y 

Z 

ACRID 

2 

3 

1.0 

2.0 

1.0 

Fi^ld; 

ID  Aerodynamic  grid  identification  number  (integer  >  0) 

CORD  Identification  number  of  ACORD  bulk  data  entry  (integer  >  0) 

Y,  Z  Location  of  the  grid  point  in  the  local  coordinate  system  defined  in 

the  ACORD  bulk  data  entry 

•  Input  Data  Entry:  AQUAD 

Description:  Defines  an  aerodynamic  surface  panel  that  has  four  corner  points 

(quadrilateral  panel) . 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

AQUAD 

ID 

IP 

IDl 

ID2 

ID3 

ID4 

AQUAD 

3 

1 

2 

3 

6 

7 
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Eield; 

ID 

IP 


Identification  number  of  the  surface  panel  (integer  >  0) 
Identification  number  of  PAERO  bulk  data  entry  for  defining  the 
aerodynamic  boundary  condition  on  the  panel  (integer  >  0) 

IDl,  ID2,  IDS,  ID4  Four  ACRID  ID*s  specifying  the  comer  points  of  the  panel 


Remarksi.. 

The  sequence  of  the  four  grid  ID*s  defines  the  out- 
normal  of  the  AQUAD  panel.  For  an  observer  standing 
on  the  wet  surface  (the  surface  on  which  the  flow  is 
passing)  of  the  panel,  the  sequence  of  these  four 
grids  must  follow  the  right-hand-rule  as  shown  in 
the  figure  to  the  right. 


IDl  ID4 


•  Input  Data  Entry:  ATRIA 

Description:  Defines  a  aerodynamic  surface  panel  that  has  three  comer  points  (triangular 

panel) . 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ATRIA 

ID 

IP 

IDl 

ID2 

IDS 

ATRIA 

4 

1 

S 

4 

7 

ID  Identification  number  of  the  surface  panel  (integer  >  0) 

IP  Identification  number  of  PAERO  bulk  data  entry  for  defining  the 

aerodynamic  bo;mdary  condition  on  the  panel  (integer  >  0) 

IDl,  ID2,  IDS  Three  ACRID  ID*s  specifying  the  corner  ^ints  of  the  panel 

Remarks: 

The  sequence  of  the  three  grid  ID*s  defines  the  out -normal  of  the  ATRIA  panel.  For  an 
observer  standing  on  the  wet  surface  (the  surface  on  which  the  flow  is  passing)  of  the  panel, 
the  sequence  of  these  three  grids  must  follow  the  right -hand-rule . 

•  Input  Data  Entzy:  AQDMEM 

Description:  Defines  a  wake  panel  on  an  arbitrary  wake  surface  that  is  shed  from  the 

trailing  edge  of  a  wing  or  the  end  section  of  a  truncated  body. 


Format  and  l^ample 


1 

2 

S 

4 

5 

6 

7 

8 

9 

10 

AQDMEM 

ID 

IP 

IDl 

ID2 

IDS 

ID4 

AQDMEM 

7 

— 

S 

4 

6 

8 

EisldJL 

ID 

IP 

IDl,  ID2,  IDS,  ID4 


Identification  nirriber  of  AQDMEM  (integer  >  0) 

Not  Used 

Four  ACRID  ID*s  specifying  the  comer  points  of  the  panel 
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Remarks: 

1-  The  figure  to  the  right  shows  an  arbitrary  wake 
surface  shed  from  the  trailing  edge  of  the  wing 
and  end  section  of  the  body.  AQDMEM  cards  are 
ertployed  to  model  this  wake  surface. 

2.  The  AQDMEM  immediately  behind  the  trailing  edge 
must  be  attached  to  the  wing/body  trailing  edge 
(i.e.  the  grid  ID  of  the  leading  edge  of  AQDMEM 
must  be  identical  to  the  grid  ID  of  the  trailing 
edge  of  the  wing/body) - 

3.  On  the  wake  surface,  the  boundary  condition 


d(l> 


=  0  is  inposed. 


dx 


•  Input  Data  Entry:  ABAR 

Description:  Defines  a  line  segment  from  which  the  wake  is  shed  downstream  to  infinity. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ABAR 

ID 

IP 

IDl 

ID2 

ABAR 

3 

— 

2 

4 

Fi^Id; 

ID  Identification  number  of  the  ABAR  bulk  data  entry 

IP  Not  Used 

IDl,  ID2,  IDS  Two  ACRID  ID*s  defining  the  line  segment 

Remarks : 

1.  If  the  wake  surface  is  assumed  to  be  a  flat  sheet  extending  from  the  trailing  edge  of  the 
wing/body,  one  ABAR  is  sufficient  to  define  one  strip  of  the  wake  sheet. 

2.  ABAR  must  be  attached  to  one,  but  no  more  than  one,  AQUAD  or  ATRIA,  (i.e.  the  line 
segment  must  be  one  of  the  sides  of  a  AQUAD  or  ATRIA  panel - 

3.  The  wake  sheet  shed  downstream  from  the  line 
segment  is  demonstrated  in  the  figure  to  the 
right.  The  edge  of  the  wake  sheet  starts  from 
the  grid  point  and  stretches  to  infinity, 
parallel  to  the  X-axis. 


4.  For  a  truncated  end  body,  the  ABAR's  must  be 
attached  to  the  trailing  edge  of  all  panels  at 
the  end  of  the  body. 


5.  In  modeling  a  thick  wing  type  of  body,  two  grid 
points  with  the  same  X,  Y  and  Z  locations  must 
be  specified  at  the  trailing  edge  of  the  body. 
Two  ABAR  *  s  are  attached  to  the  upper  and  lower 
side  of  the  trailing  edge.  In  this  way,  the 
potential  jump  of  the  wake  effect  can  be 
represented  by  the  potential  difference  between 
these  two  wake  sheets. 
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*  Input  Data  Entzy:  PAEP0 

Description:  Defines  the  aerocfynaniic  boundary  condition  of  a  surface  panel. 


Format  and  Example 


1 

3 

4 

5 

6 

7 

8 

9 

10 

PAERO 

ID 

FLWRT 

PAERO 

3 

0.0 

Eisldl 

ID  Identification  number  of  PAERO 

FLWRT  Amount  of  mass  ratio  (in  percentage)  which  is  allowed  to  penetrate 

into  the  panel 

FLWRT=0.0  represents  a  solid  surface  (the  zero-normal-velocity 
boundary  condition  is  inposed) 

FLWRT=100.0  represents  an  inlet  face  with  no  flow  leakage 

•  Input  Data  Entzy:  AOOED 

Description:  Defines  a  group  of  aerodynamic  conponents. 


Format  and  Example 


1 

3 

4 

5 

6 

7 

8 

9 

10 

ACORD 

ID 

CP 

DELTA 

XMCNT 

YMCNT 

ZMCNT 

ACORD 

3 

7 

0.0 

1.0 

2.3 

4.0 

Field: 

ID 

CE> 

DELTA 


XMCNT,  YMCNT,  ZNCNT 


Identification  number  of  the  ACORD  bulk  data  entry  (All  AGRID  data 
entries  must  refer  to  ACORD) 

Identification  number  of  local  coordinate  system  in  which  the 
location  of  AGRID  is  defined 

Inclination  angle  of  the  body  axis  to  the  free  stream  defined  in  the 
local  coordinate  system,  positive  nose  up,  and  measured  in  degrees. 
This  parameter  will  not  change  the  panel  position.  Its  effect  is 
introduced  in  the  boundary  condition.  Therefore,  DELTA  must  be  a 
small  value.  Since  most  under-wing  stores  have  a  small  inclination 
angle  to  the  free  stream,  this  input  would  facilitate  the  user  with  a 
sinpler  definition  of  inclined  bodies 

Pitch  and  yaw  moment  center  defined  in  the  local  coordinate  system 
and  used  only  for  calculating  the  pitch  and  yaw  moment  of  each 
subsystem 


Lifting  Surface  Input 


In  defining  an  aerodynamic  panel  element  (a  lifting  surface),  the  CAEROl  data  entry  of 
ASTROS  is  modified  to  cope  with  the  thickness  distribution  for  Z0NA51U  and  steady  pressure 
distribution  for  ZTAIC.  The  connectivity  of  the  lifting  surface  to  a  surface  panel  model  can  also 
be  specified.  This  is  done  by  defining  a  new  bulk  data  entry  "CAER07"  as  follows: 


•  Input  Data  Entzy:  CAERD7 

Description:  Defines  a  lifting  surface  (panel)  for  the  Z0NA6,  Z0NA7,  Z0NA51U  and  ZTAIC 

methods. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CAER07 

EID 

IP 

NSPAN 

NCHORD 

LSPAN 

ISCP 

+ABC 

CAER07 

1 

1 

10 

10 

3 

2 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

C3^07 

XRL 

YRL 

ZRL 

RCH 

LRCHD 

RFOIL 

+AB 

+ABC 

0.0 

0.0 

0.0 

1.0 

3 

1 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

CAER07 

XTL 

YTL 

ZTL 

TCH 

LTCHD 

TFOIL 

+AB 

0.0 

1.0 

0.0 

0.5 

10 

1 

Field; 

EID 

IP 

NSPAN 

LSPAN 

ISCP 

XRL,  YRL,  ZRL 

RCH 

LRCHD 


RFOIL 

XTL,  YTL,  ZTL 

TCH 

LTCHD 

TFOIL 


Identification  nxairiber  of  CAER07 

Identification  number  of  ACORD  entry  specifying  the  group  of  the 
CAER07 

Number  of  spanwise  boxes 
Number  of  chordwise  boxes 

Identification  number  of  ITAIC  entry  for  specifying  the  steady 
pressure  input  of  each  strip  (Used  only  for  the  ZTAIC  method) 

Location  of  the  leading  edge  at  root  in  the  local  coordinates  and 
defined  in  the  ACORD  entry 
Root  chord  length 

If  LRCHD>0,  LRCHP  is  the  identification  number  of  an  AEFACT  data 
entry  containing  a  list  of  division  points  for  chordwise  boxes  along 
the  root  chord 

If  LRCHD<0,  -LRCHD  is  the  identification  number  of  an  AEFACT  data 
entry  containing  a  list  of  ACRID  ID*s.  Thus  the  root  is  attached  to 
a  surface  panel  model.  This  provides  the  "connectivity"  information 
of  the  lifting  surface  model  with  the  surface  panel  model. 
Identification  of  a  PAFOIL  entry  specifying  the  thickness  and  camber 
of  the  airfoil  at  the  root  (Only  used  for  the  Z0NA51U  method) 

Same  as  XRL,  YRL,  ZRL,  but  for  the  tip 

Tip  chord  length 

Same  as  LRCHD,  but  for  the  tip 

Same  as  RFOIL,  but  for  the  tip 


•  Input  Data  Entry:  ITAIC 

Description:  Defines  a  set  of  AEFACT  entries  for  steady  pressure  input  for  each  strip. 


Format  and  Example 


1 

3 

4 

5 

6 

7 

8 

9 

10 

ITAIC 

ISCP 

NXI 

lAXl 

lACPUl 

lACPLl 

NX2 

IAX2 

IACP02 

+ABC 

ITAIC 

4 

20 

1 

3 

3 

20 

1 

4 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

ITAIC 

IACPL2 

— 

— 

— 

NXN 

lAXN 

lACPUN 

lACPLN 

+ABC 

+ABC 

3 

— 

— 

— 

21 

10 

1 

2 

ISCP 

NXi 

lAXi 

lACPUi 

lACPLi 


Identification  number  of  ITAIC  indexed  by  the  CAER07  entry  if  the 
ZTAIC  method  is  used 

The  number  of  steady  pressure  points  to  be  input  along  the  i-th  strip 
ID  of  AEFACT  entry  that  defines  NXi  number  of  X  locations  of  the 
steady  pressure  (0.0<X<1.0)  of  the  i-th  strip 

ID  of  AEFACT  entry  that  defines  NXi  nvimber  of  pressure  coefficients 
on  the  upper  surface  of  the  i-th  strip 

ID  of  AEFACT  entry  that  defines  NXi  number  of  pressure  coefficients 
on  the  lower  surface  of  the  i-th  strip 


The  total  nxmiber  of  NXi,  lAXi,  lACPUi  and  lACPLi  must  be  NSPAN,  where  NSPAN  is  defined  in  the 
CAER07  entry. 
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Input  Data  Entzy:  PAFOIIi 


Description:  Defines  the  airfoil  thickness  and  camber  for  the  Z0NA51U  method. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

PAFOIL 

PID 

NX 

lAFX 

lATU 

lATL 

ICAM 

RADIUS 

PAFOIL 

PID 

21 

3 

4 

7 

5 

0.10 

Eieldj. 

PID  Identification  number  of  PAFOIL  indexed  by  the  GAER07  entry 

NX  Number  of  X  locations  in  percent  chord  to  define  the  thickness  and 

camber  (0.0  <  X  <  100.0) 

lAFX  Identification  number  of  AEFACT  entries  containing  NX  numbers  of  X 

locations.  The  first  and  last  values  must  be  0.0  and  100.0, 
respectively . 

lATU  Identification  number  of  AEFACT  entries  containing  NX  numbers  of 

upper  surface  ordinates  in  percent  chord 
lATL  Same  as  lATU,  but  for  the  lower  surface 

ICAM  Identification  number  of  AEFACT  entries  containing  NX  numbers  of 

camber  ordinates  in  percent  chord 
RADIUS  Radius  of  leading  edge  in  percent  chord 

It  is  sufficient  that  all  geometric  data  needed  by  the  ZONA6,  ZONA7,  ZONA51U  and  ZTAIC 
codes  can  be  provided  by  the  new  bulk  data  entries.  We  believe  that  the  proposed  AGM  is  also 
capable  of  coping  with  other  higher-order  panel  codes  if  they  are  incorporated  into  ASTROS. 

3.2  Generation  Of  Aerodynamic  Influence  Coefficient  (AIO  Matrices 

The  primary  function  of  the  ZAERO  module  is  to  generate  Unified  Aerodynamic  Influence 
Coefficient  (UAIC)  matrices  for  the  entire  Mach  range.  Since,  by  definition,  the  AIC  matrix  is 
independent  of  the  structural  deformations  (i.e.  the  modes),  the  ZAERO  module  will  be 
incorporated  outside  of  the  ASTROS  optimization  loop.  All  AIC  matrices  are  computed  only 
once  and  saved  in  the  CADDB  database. 

Once  all  of  the  AIC  matrices  are  generated  by  2^AERO,  the  Generalized  Aerodynamic  Forces 
(GAFs)  for  each  optimization  step  have  to  be  computed.  Since  the  UAICs  GAF's  share  the 
same  definition  as  those  generated  by  the  Doublet  Lattice  Method  (DLM)  and  Constant  Pressure 
Method  (CPM),  the  GAFs  can  be  adopted  directly  by  ASTROS’  p-k  method  for  flutter 
analysis/design. 

In  what  follows,  the  AIC  matrices  of  the  ZAERO  module  (including  ZONA6,  ZONA7,  ZTAIC 
and  ZONA51U)  are  formulated.  In  the  absence  of  a  surface  panel  model  these  AIC  matrices  are 
equivalent  the  "QKK"  matrices  of  ASTROS  generated  by  the  DLM  or  CPM  lifting  surface 
methods.  However,  for  surface  panel  methods,  diese  AIC  formulations  are  generalized  for  wing- 
body  configurations  and  are  discussed  as  follows: 

Influence  Coefficient  Matrices  For  Wing-Body  Configurations 

As  discussed  in  the  Aerodynamic  Geometry  Module  (AGM),  ZONA6  and  ZONA7  are  the  hybrid 
methods  in  which  the  aircraft  configuration  is  broadly  divid^  into  two  categories:  the  body-like 
components  modeled  by  surface  panel  methods  and  the  wing-like  components  by  lifting  surface 
methods.  A  constant  unsteady  source  singularity  (a)  is  applied  on  each  body  panel,  while,  a 
constant  unsteady  pressure  singularity  (ACp)  is  distributed  on  each  wing  panel.  Five  influence 
coefficient  matrices,  namely  [PIC],  [UIC],  [VIC],  [WIC]  and  [NIC]  for  the  potential  axial 
velocity,  lateral  velocity,  vertical  velocity  and  normal  velocity  influence  coefficient  matrices. 
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respectively,  are  constructed  to  relate  the  unknown  singularities  (cr  and  ACp )  to  the  perturbation 
potential  and  velocities  {u},  {v}  and  {w}  at  all  panels,  i.e. 


I4  =  [«C]| 

Ucp 


(«)=[C//C](  )  = 

\acpI 

{v)=(WC](  ]-- 

Ucpj 

[w]=[WIC][  °  ] 
yACp) 


PICbb  PICwb  I  o 
PICew  PICwwl  {ACp 


UICbb  UICwb]  a  I 
UICbw  UICww\\^Cp\ 

VICbb  VICwb  1  /  o-  \ 
VICbw 

WICbb  WICwb]!  a 
[wICbw  WICwwWac 


(3.2.1) 


(3.2.2) 


(3.2.3) 


(3.2.4) 


where  n  =  n^i  +  /y  +  rizk  is  the  out-normal  vector  of  each  panel  and  the  subscripts  denote: 

BB  =  the  influence  at  the  body  control  points  due  to  the  body  panels 

BW  =  the  influence  at  the  wing  control  points  due  to  the  body  panels 

WB  =  the  influence  at  the  body  control  points  due  to  the  wing  panels 

WW  =  the  influence  at  the  wing  control  points  due  to  the  wing  panels 

Notice  that  unlike  the  lifting  surface  methods,  such  as  the  Doublet  Lattice  Method  (DLM)  and 
Constant  Pressure  Method  (CPM)  in  which  only  one  influence  matrix  is  required,  namely,  the 
partition  matrix  NIC^^  in  ^  (3.2.5),  the  wing-body  configurations  require  five  matrices  for  the 
unsteady  force  computations.  The  unknown  singularities  cr  and  ACp  are  solved  by  relating  the 
right  hand  side  of  Eq  (3.2.5)  to  the  wing-body  boundary  conditions  which  will  be  discussed  next. 

Unsteady  Wing-Body  Boundary  Condition 

The  boundary  condition  of  the  wing-tike  component  (lifting  surface  method)  is  the  well  known 
expression: 


(m,  V,  w )  •  n  =  Fw 


(3.2.6) 


where  Fw  =  —  -i-  ikh 
dx 

h  is  the  normal  displacement  (or  the  mode  shape) 

Bh 

—  is  the  slope  of  h 
dx 

Eq  (3.2.6)  is  the  standard  boundary  condition  employed  for  the  lifting  surface  methods  such  as 
DLM,  CPM,  ZONA51U  and  the  ZTAIC  codes.  For  an  arbitrary  body-like  component 
experiencing  oscillatory  motion,  the  body  fixed  boundary  condition  on  ±e  body  surface  reads: 
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[v  -  VB)*n  =  0  at  So  =  0 


(3.2.7) 


where 


^  =  So(  X,  y,  z  ),  the  mean  body  surface  function 

V  =  total  flow  velocity  on  the  body  surface  So  =  0 
Vb  =  the  velocity  due  to  the  body  surface  motion 

The  elastic  body  modes  can  be  generalized  to  contain  two  sets  of  components,  namely,  the  one 
defined  in  the  pitch  plane  h^,  and  the  one  defined  in  the  yaw  plane  hy.  Perturbing  Eq  (3.2.7)  with 
respect  to  the  small  elastic  body  modes  and  hy  results  in  a  generalized  unsteady  boundary 
condition  in  the  body-fixed  system: 


{u,v,w)»n  =  FB 

(3.2.8) 

where 

Fb  =  -ZBfixihzUo  +  ikhz)  -  rizik 

z  +  ik  hz) 

-yBTixihyUo  +  ikhy)  -  ny[h 

■Y  +  ik  hr] 

(3.2.9) 

2b  ,  JB 
Uo. 

hz,  hr 
hz,  hy 

tt  II 

hz  ,  hy 
( )’.  ( )" 


the  z  and  y  coordinate  of  the  control  point 

the  steady  perturbation  velocity  in  the  jc-direction 

body  mode  shapes  projected  in  the  x-z  andx-y  plane  respectively 

the  slopes  of  hz  and  hy  respectively 


the  curvatures  of  hz  and  hy  respectively 
d  dP" 

~  ^  (  )»  differential  operator  with  respect  to  x. 


Note  that  in  the  boundary  condition  for  lifting  surfaces,  Eq  (3.2.6),  the  downwash  function  Fw 
is  only  a  function  of  the  wing  mode  h  and  the  reduced  fi-equency  k.  By  contrast,  the  downwash 
function  Fb.  of  Eq  (3.2.9)  for  bodies  is  a  function  of  normal  vector  n  ,  the  steady  velocity,  the 
mode  shapes  hy  and  hz  and  their  derivatives,  and  reduced  frequency  k.  This  implies  that  while 
the  unsteady  lifting  surface  solution  is  totally  uncoupled  from  the  steady  mean  flow  influence, 
the  unsteady  body  solution  must  include  the  steady  mean  flow  influence  which  enters  through 
the  boundary  condition,  Eq  (3.2.9). 

In  contrast  to  the  Fw  expression  in  Eq  (3.2.6),  the  F^  expression  in  Eq  (3.2.9)  contains  second 
order  derivatives  in  hz  and  hy,  which  are  related  through  a  required  transformation  from  the  wing- 
fixed  coordinate  system  to  the  present  body-fixed  coordinate  system. 

Combining  Eqs  (3.2.5),  (3.2.6)  and  (3.2.8)  yields  the  unknown  singularity  strength  on  the  wing¬ 
like  and  body-like  components,  i.e. 


<J 

ACp 


(3.2.10) 
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Pressure  Coefficients,  Generalized  Forces 

Based  on  the  work  of  Garcia-Fogeda  and  Liu  (Ref  45),  the  unsteady  pressure  coefficient  on  the 
body-like  component  expressed  in  the  body-fixed  coordinates,  reads 


Cp 


where 


So 


-2So{{  1  +Mo)( «  -  h'z  Uo  -  yB  Jty  Uo  )  +  v^v  +  WoW 

+  h^Wo  +  hyVo  +  ik{(l>  +  hzWo  +  hyVo  -  zskluo  -  yshyUo] 


1- 


y-1 


mI(2uo 


+vi  + 


y-l 


(3.2.11) 


Uo,  Vo,  Wo  =  the  steady  mean  perturbation  velocities  on  the  body 

^,u,v,w  -  the  unsteady  perturbation  potential  and  velocities  on  the  body 


According  to  Eq  (3.2.11),  the  unsteady  pressure  on  body  involves  the  steady  solution  of 
Uo,  Vo  and  Wq.  In  principle,  Eq  (3.2.11)  is  also  applicable  to  lifting  stufaces.  When  placed  at 
zero  degrees  angle  of  attack,  the  steady  velocities  for  lifting  surfaces  are  =  w©  =  0. 

Eq  (3.2. 1 1)  is  then  reaiuced  to  the  well-known  expression. 


ACp  -  -2  {u  +  ik  (3.2.12) 

It  should  be  remarked  that  the  unsteady  pressure  Cp  for  bodies,  Eq  (3.2.1 1),  involves  coupling 
terms  with  the  perturbation  velocities  of  the  steady  mean  flow,  which  incorporates  the  body 
thickness  effect.  By  contrast,  the  unsteady  pressure  ACp  for  lifting  surface,  Eq  (3.2.12)  is 
imcoupled  from  the  steady  mean  flow  term. 

The  generalized  aerodynamic  forces  (GAFs),  defined  as  the  work  done  by  the  unsteady  forces, 
can  be  expressed  in  terms  of  the  mode  shapes  and  the  pressure  coefficients,  i.e. 


[/  \  /  \ 


i  =  1 
NW 

+  ^  ACpf^Aw,hl 
»■  =  1 


(/) 


(3.2.13) 


where 


cj,r 

=:  1 

AC^pP 

=  1 

Ab, 

=  1 

Aw, 

=  *1 

hiW  h 

U)  _ 

Yi  - 

dhPP 

dx  ’ 

dx 

th&i-th  body  panel 
=  the  I-th  mode  shape  on  the  i-th  wing  panel 
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NB,  NW  are  the  number  of  body  panels  and  wing  panels,  respectively 

When  the  generalized  forces  for  the  body  component  (first  summation  terms  of  Eq  (3.2.13) )  are 
compared  with  that  for  the  wing  component  (the  second  summation  term),  it  can  be  seen  that  the 
terms  associated  with  the  I-th  mode  in  Qjj  due  to  the  body  thickness  effect  are  essentially 


5 


dhP'>  dhP^ 


‘dx 


-ZBi  +  nxi 


dx 


This  term  is  not  accounted  for  in  the  slender-body  approach  employed  in  the  Doublet  Lattice 
Method. 


Solution  Procedure  For  Computing  Unsteady  Pressure  And  Generalized  Forces  Of  Wing-Body 
Configurations 

Fig  (3.2.1)  presents  the  solution  procedure  in  obtaining  the  wing-body  unsteady  pressure  and 
generalized  forces,  which  develops  as  follows: 


Figure  3.2. 1  Flow  Chart  of  ZONA6/ZONA7  Computation  Procedure. 
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a)  Solve  for  steady  perturbation  velocities  Uo,  Vq  and  Wo. 

b)  Construct  Influence  Coefficient  matrices  [  i7/C  ],  [  VIC  ],  [  WIC  ],  [  PIC  ]  and[  NIC  ]. 

c)  Construct  Downwash  functions  Fb  and  Fw  of  Eqs  (3.2.6)  and  (3.2.9);  perform  matrix 
decomposition  on  [  iV/C  ]  and  solve  for  <t  znd.  ACp  from  Eq  (3.2.10). 


d)  Compute  unsteady  velocities  u,  v  and  w  and  unsteady  potential  (j>  according  to  Eqs  (3.2.1) 
through  (3.2.4). 

e)  Compute  unsteady  pressure  Cp  on  body  according  to  Eq  (3.2. 1 1). 

f)  Compute  generalized  forces  Qj j  according  to  Eq.  (3.2. 13). 

The  above  procedures  require  the  structural  modes  in  the  computations  of  the  downwash 
functions  Fg  andF^^^,  the  unsteady  pressure  computation  of  Cp  and  generalized  forces  (2//-  This 
implies  that  these  computations  must  be  within  the  design  loop  of  ASTROS  which  would  be 
computationally  inefficient.  Therefore,  reformulation  of  Ae  above  procedures  to  construct  a  set 
of  AIC  matrices  which  are  independent  of  the  structure  is  required  for  ASTROS 
design/optimization. 

Construction  Of  Unified  Aerodynamic  Influence  Coefficient  (UAIC)  Matrices 

The  definition  of  AIC  is  a  square  matrix  which  relates  the  downwash  to  the  unsteady  and  lifting 
forces.  To  constract  an  AIC  matrix  for  wing-like  components  is  rather  straightforward.  It  can  be 
obtained  according  to  the  following  derivation: 


=  [[0]  .  [A»J][W/C  ]•' 


(3.2.14) 


where 


=  the  normal  forces  on  wing  panels 


[  Aw  J 


=  a  diagonal  matrix  in  which  the  element  Aw^  represents  the  area  of  the 
i-th  wing  panel 


The  procedure  to  obtain  AIC  for  the  body-like  components  is  more  complicated  than  that  for 
wing-like  components.  Rewriting  Eq  (3.2.11)  in  matrix  form  and  multiplying  by  the  body  panel 
area  yields: 


{Lbody  ) 


[  Bbb  ]  ,  [  Bwb  ] 


where 


[  [Rbb]  ,  [Bwb]  ] 


r  ab  \\-2So\ 


[  1+MoJ 


[  UICbb  ] 


[  UICWB  ] 


(3.2.15) 


+  [v.J 


[  VICbb  ] 


[  VICwB  ] 
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+  [WoJ 


[  WICbb  ]  ,  [  WICwB  ] 


+  [  J 


[PICbb]  ,  [PICwb] 


[MC]-1 


(3.2.16) 


{L-body }  =  the  normal  forces  on  body  panels, 
d  =  -IAS o[-{l+Uo)(  zbI^Uo  +  yB  hyUo) 

+  hzWo  +  hy  Vo  +  i^(  hzM>o  +  hyVo  -  ZBh'zUo  -  yBhyUo) 

(3.2.17) 

I  Ab  J  =  the  element  represents  the  area  of  the  i-th  body  panel 

Combining  Eqs  (3.2.14)  and  (3.2.15)  and  separating  the  modes  and  the  influence  coefficients 
reveals  that  six  matrices  are  required  to  relate  the  modes  to  the  panel  forces. 


[Lbody 

J^ing 


=  [  A/Cz  ] 


f  hz 


+  [  A/Czx  ] 


+  [  AICzxx  ] 


d^hz 

dx'^ 

0 


+  [  A/Cr  ]  +  [  AlCyx  ] 

[Oj 


dx 


1  ® 


+  [  A/Cm  ; 


d^hr 

dx'^ 

0 


(3.2.18) 


where  [A/C^],  [AICzx\  and  [A/CzxjJ  are  the  AIC  matrices  that  relate  the  vertical  displacements 
and  their  slopes  and  curvatures  to  the  forces  on  the  body  and  wing  components  respectively, 
whereas  [A/Cy],  [AICyx]  and  [AICyxx  ]  are  for  the  lateral  motion. 


The  construction  of  the  six  AIC  matrices  requires  a  major  modification  of  the  current 
ZONA51U,  ZONA6,  ZONA7,  ZTAIC  codes.  The  ASTROS  matrix  utility  routines  will  be  used 
for  the  matrix  manipulations  of  Eqs  (3.2.14)  through  (3.2.18). 

Spline  Matrix  Of  Wing-Body  Configurations 

Eq  (3.2.18)  shows  that  in  addition  to  the  modes  and  slopes,  the  curvatures  are  also  required  in 
computing  the  wing-body  unsteady  aerodynamic  forces.  The  surface  spline  technique  in 
ASTROS  can  be  adapted  directly  for  the  modal  interpolations  of  lifting  surface  models  (i.e. 
wing-like  components).  Similarly,  the  beam  splme  technique  of  ASTROS  can  be  adapted 
directly  for  the  interpolations  of  modes  and  slopes  on  surface  panel  models  (i.e.  body-like 
components).  However,  the  interpolation  technique  of  curvatures  is  lacking  in  the  current 
version  of  ASTROS.  A  new  beam  spline  technique  will  be  developed  that  adopts  the  current 
ASTROS'  beam  spline  technique  as  a  basis  and  computes  the  curvature  through  analytical 
differentiation  of  the  slopes. 

Once  the  new  beam  spline  technique  is  developed,  spline  matrices  can  be  defined  to  relate  the 
modal  data  at  the  structural  grids  to  the  modes,  slopes  and  curvatures  at  the  aerodynamic  grids. 
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Consequently,  the  six  AIC  matrices  defined  in  Eq  (3.2.18)  can  be  transformed  into  the  stmc^ral 
grids  by  substituting  these  spline  matrices  into  equation  (3.2.18)  and  the  generalized 
aerodynamic  forces  can  be  computed  based  on  equation  (3.2.13)  accordingly. 

ASTR0SI2AER0  Integration 

Under  the  proposed  effort,  four  new  methods  to  generate  the  unsteady  aerodyn£umc  influence 
coefficients  for  use  in  ASTROS’  ASE  analyses  will  be  implemented.  The  original  Doublet 
Lattice  Method  (DLM)  and  Constant  Pressure  Method  (CPM)  will  be  modified  to  conform  to  the 
new  aerodynamic  matrix  formats  which  will  be  required  to  support  the  archived  import  features 
that  are  proposed  for  all  AIC  generation  methods. 

For  integration  into  ASTROS,  two  complexities  have  been  introduced  that  require  modifications 
to  ASTROS.  First,  with  the  advent  of  numerous  methods  in  each  Mach  range,  the  user  must  tell 
ASTROS  which  method  to  use.  In  the  past,  the  Mach  number  essentially  played  that  role. 
Secondly,  ASTROS  must  identify  the  AIC's  by  more  than  the  Mach  number,  reduced  frequency 
and  symmetry  (M,  k,  sym)  triplet  that  is  used  for  DLM  and  CPM.  Now  the  method  (ZONA51U, 
ZONA6,  ZONA7  and  ZTAIC)  and  splining  strategy  must  be  stored  as  well.  In  order  to  make  the 
grouping  more  straightforward,  and  to  allow  the  ARCHTVE/IMPORT  features  to  support  ASE 
response  analyses,  ASTROS  will  be  modified  to  allow  the  user  to  name  a  collection  of  input 
describing  a  discretized  aerodynamic  configuration.  That  named  model  can  then  be  ARCHIVED 
or  IMPORTED  as  a  single  entity  (i.e.  a  single  collection  of  database  entities). 

To  support  the  new  methods  in  a  more  straightforward  manner,  the  proposed  AGM 
(Aerodynamic  Geometry  Module)  will  be  implemented  to  create  a  single  form  (set  of  database 
entities)  that  is  capable  of  describing  all  the  various  parts  of  any  of  the  supported  methods. 
Those  geometrical  descriptions  that  are  common  will  be  indistinguishable,  with  additions  and 
modifications  handled  using  rationalized  database  stractures  (relations  and/or  matrices). 

The  existing  DLM  and  CPM  in  ASTROS  will  be  modified  to  require  their  geometrical 
description  in  the  same  form  as  that  of  the  AGM.  This  implies  that  the  UNSTEADY  module  in 
ASTROS  will  be  broken  apart  into  a  geometry  engine  and  a  computational  engine.  These 
improvements  would  become  permanent  features  in  standard  ASTROS. 

Another  ramification  of  the  ZAERO  module  is  that  the  splining  methodology  will  become 
method-dependent.  Lifting  surfaces  in  ASTROS  are  already  treated  differently  than  bodies  -  the 
surface  spline  is  only  available  for  lifting  surfaces.  While  the  basic  spline  technologies  (beam 
spline  and  infinite  stuface  spline)  are  adequate,  their  application  will  become  method  dependent. 
This  implies  that  the  ASTROS  low  level  spline  code  will  be  used  within  an  augmented  SPLINE 
module  that  will  handle  the  new  requirements.  First  among  the  necessary  chmges  is  to  make  use 
of  the  AGM  geometrical  descriptions  rather  than  the  current  ASTROS  capabilities. 

As  a  potential  follow-on  effort,  the  AGM  geometiy  (which  is  very  similar  to  an  FE  mesh)  would 
be  made  available  in  forms  compatible  with  graphical  pre/post  processors  to  display  aerodynamic 
output  quantities  graphically.  These  features  Income  readily  available  once  the  AGM  and  AIC 
data  are  archived  to  a  CADDB  database.  External  software  could  be  developed  by  the  end-user 
or  third  parties  to  support  AIC  model  development  in  a  very  general  way.  This  potential  offers  a 
chance  to  resolve  one  of  the  major  limitations  of  ASE  an^ysis  in  commercial  FE  codes  in  that 
the  FE  codes  do  not  support  "model  checkout"  and  debug  features  for  the  unsteady  aerodynamic 
model. 
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SECTION  4 


SOFTWARE  DESIGN  BLUE-PRINT  OF  ASTROS/ASE  MODULE 
4.1  Scenarios  for  ASTROS/ASE  Applications 

As  discussed  previously  in  section  D,  the  proposed  ASTROS/ASE  module,  once  integrated  with 
ASTROS,  will  bring  to  the  ASTROS  MDO  environment  a  multidisciplinary  synergistic 
technology  that  will  provide  a  quick  turnaround  aeroelastic  design/analysis  capability  tiuroughout 
the  complete  flight  regime.  The  following  scenarios  together  with  full-scale  case  studies  of  the 
AAW  will  elucidate  the  more  versatile  features  of  the  ASTROS/ASE  capability. 

Active  Flexible  Wing  (AFW)  Roll  Performance  with  Flutter  Constraints 


Figure  4.1.1 

Figure  4.4.1  presents  an  AFW  wind  tunnel  model  which  has  two  leading  edge  and  two  trailing 
edge  control  surfaces  to  meet  the  roll  performance  requirements.  An  initial  control  system 
design  satisfies  these  constraints  but  introduces  antisymmetric  flutter.  With  the  proposed 
ASTROS/ASE  module,  the  control  gains  of  the  four  control  surfaces  are  added  to  the  list  of  the 
design  variables.  The  new  ASE  module  and  the  static  aeroelastic  module  are  employed  to 
optimize  the  control  gains  and  the  structural  design  variables  for  satisfying  the  flutter  and  roll 
performance  constraints  with  a  minimal  weight  penalty. 

Trade-Off  Study  of  Gust  Load  Alleviation  by  Passive  and  Active  Means 
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A  drone  aircraft  such  as  the  DAST  (Drones  for  Aerodynamic  and  Structural  Testing,  Ref  46) 
shown  in  figure  4.1.2,  is  experiencing  excessive  wing  root  RMS  bending  moment  in  open-loop 
continuous  gust  analysis.  With  the  new  ASE  module,  the  baseline  structure  is  used  as  a  starting 
point  for  two  design  studies: 

•  Passive  Design  -  Add  gust  response  constraints  and  repeat  the  structural  design  process. 

•  Active  Design  -  Extract  the  plant  state-space  model  for  designing  a  load  alleviation 

control  system  outside  of  ASTROS. 

-  Due  to  the  control  system  limitations,  the  control  law  helps  but  does  not 
fully  satisfy  the  gust  response  requirements. 

-  Redesign  the  structure  in  the  presence  of  the  control  system  with  gain- 
margin,  phase  margin  and  singular  value  constraints  added  to  the  list  of 
constraints. 

The  results  of  both  of  the  design  are  compared  for  trade-off  studies  of  structural  changes  versus 
complexity  of  the  control  system. 

Store-Flutter  Suppression  of  a  Generic  Advanced  Fighter  (GAF)  with  Structural 
Uncertainties 


Figure  4.1.3 

Figure  4.1.3  shows  a  generic  advanced  fighter  (GAF)  with  various  store  carriages  for  the  air-to- 
ground  mission,  requiring  flutter  clearance  assurance.  However,  uncertainties  of  the 
stores/pylons  structural  properties  may  introduce  an  unexpected  flutter  problem  which  could  be 
resolv^  by  a  robust  control  system  for  store  flutter  suppression.  The  state-space  matrices  of  the 
aircraft  with  external  stores  and  the  derivatives  of  these  matrices  with  respect  to  the  structural 
properties  of  the  stores/pylons  will  first  be  generated  by  the  new  ASE  module.  MATLAB's  p.- 
analysis  tool  box  will  then  be  used  to  design  a  control  system  which  suppresses  flutter  over  a 
range  of  parameter  uncertainties.  The  designed  control  law  will  be  introduced  to  ASTROS  for 
the  verification  with  various  store  models  at  different  flight  conditions. 

Full-Scale  Active  Aeroelastic  Wing  (AAW)  Studies 

The  concepts  of  Active  Aeroelastic  Wing  (AAW)  Technology  envision  the  use  of  wing 
aeroelastic  flexibility  as  an  advantage,  with  the  wing  control  surfaces  employed  as  tabs  to 
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wing  twist  for  aer^^amic  control,  whereas  the  current  practice  is  to  increase  the  wing 
stiffness  in  order  to  minimize  flexibility.  The  AAW  could  offer  a  wing  system  that  is 
si^ificantly  lighter  in  weight  than  current  supersonic  fighters  while  improving  controllability 
and  aerodynamic  efficiency.  AAW  technology  has  been  demonstrated  through  analysis  and 
wmd  tunnel  modeling  (Refs  1, 2),  and  a  full  scale  flight  demonstration  will  be  performed  in  the 
ne^  foture  as  evidenced  by  the  recent  PRDA  firom  the  Air  Force.  While  keeping  the  wing 
sufficiently  flexible,  the  Active  Aeroelastic  Wing  may  be  subjected  to  a  number  of  dynamic 
aeroelastoc  problems  such  as  flutter,  gust  response,  and/or  store  flutter.  Apparently,  flutter 
suppression  and  gust  load  alleviation  by  active  means  would  be  the  feasible  solution,  which  calls 
for  a  mu^sciplinary  design/analysis  capability  such  as  the  proposed  ASTROS/ASE  software 
system,  ^e  ZONA  team  will,  by  utilizing  the  ASTROS/ASE  module,  perform  a  parallel  study 
of  the  i^W  full  scale  model  design/analysis.  The  outcome  of  the  full  scale  flight  demonstration 
could  also  serve  as  a  validation  test  case  for  the  ASTROS/ASE  module.  ZONA  Technology  is 
currently  working  out  with  North  American  Rockwell/Seal  Beach,  CA,  a  MOU  concerning  the 
possible  data  exchange  and  collaboration  on  the  proposed  AAW  project. 

The  following  sections  layout  a  software  design  blueprint  of  the  proposed  ASTROS/ASE 
module.  A  detailed  flowchart  of  the  ASE  module  that  depicts  all  of  its  submodules  and  their 
interrelationships  is  represented.  The  functions  of  each  submodule  and  its  input  and  output  in 
terms  of  matnx  entities  are  discussed.  Inteirelationships  between  submodules  are  also  identified 
in  teims  of  ASTROS/MAPOL  sequences  and  the  CADDB  database  management  system.  New 
bulk  data  entries  are  defined  as  required. 

4.2  Overview  of  the  ASE  Module 

The  planned  ASE  module  will  facilitate  the  inclusion  of  multi-input  multi-output  (MIMO) 
conttol  system  effects  on  the  dynamic  stabihty  and  response  in  the  ASTROS  multidisciplinary 
analysis  and  design  optimization  software  package.  Its  overall  capabilities  will  include: 

(i)  Provide  closed-loop  robust  stability  analysis. 

(ii)  Add  continuous  gust  response  capabilities. 

(iii)  Allow  the  inclusion  of  stability  and  gust-response  constraints  in  structural  design 

optimization.  ^ 

(iv)  Allow  the  inclusion  of  user-defined  control  parameters  of  a  given  control  law  in  the 
multidisciplinary  optimization  process. 

(v)  Export  an  efficient  state-space  representation  of  the  aeroservoelastic  system  for 
subsequent  analysis  and  control  synthesis  with  commercially  available  tools  such  as 
MATLAB  and  MATRIX  x. 

The  ASE  module  vstill  be  based  on  state-space  formulations.  The  structure  is  represented  by  a  set 
of  baseline  normal  modes  serving  as  generalized  coordinates.  The  unsteady  aerodynamic  forces 
wiU  be  represented  by  minimum-state  rational  approximations  of  the  ZAERO  module  generated 
ttanscendental  frequency  domain  generdiz^  force  coefficient  matrices.  The  control  system  will 
be  represented  by  a  state-space  realization  of  user-defined  series  of  polynomial  transfer 
unctions.  A  gust  filter  will  be  defined  such  that  a  white-noise  input  produces  an  approximation 

of  either  Dryden’s  or  von  Karman’s  power  spectral  density  (PSD)  of  atmospheric  continuous 
gusts. 

Tlie  stability  analysis  and  constraints  will  be  based  on  root-loci,  Nyquist  curves  and  transfer- 
functton  singular  values  in  the  frequency  domain.  The  gust  response  analysis  and  sensitivities 
will  be  based  on  the  stochastic  Lyapunov  formulation. 

A  physical  weighting  algorithm  will  rate  the  aerodynamic  data  terms  according  to  their 
aeroelastic  importance.  These  rates  will  be  used  for: 


37 


STTR  Final  Rpt050S96 


(i)  Weighting  the  data  terms  in  the  rational  approximation  process. 

(ii)  Rating  the  modes  for  application  of  size-reduction  techniques. 

There  will  be  several  options  for  the  reduction  of  the  order  of  the  state-space  equations.  These 
options  will  allow  a  combination  of  modal  truncation,  static  residualization,  and  dynamic 
residualization. 

The  new  features  will  be  applicable  to  open-loop  as  well  as  closed-loop  systems. 


Figure  4.2.1  General  Flow  Chart  of  the  ASE  Module. 
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4.3  Description  Of  The  ASE  Blocks 


A  flow  chart  of  the  computation  sequence  associated  with  the  planned  ASE  capabilities  is  shown 

in  Figure  4.2.1.  The  flow  chart  is  divided  into  three  parts; 

(i)  Preface  structiu^  and  aerodynamic  blocks  which  construct  the  g-set  finite-element 
structural  model  and  the  ^-set  aerodynamic  panel  model,  and  the  spline  matrices  that 
convert  displacements  and  slopes  between  the  two  sets. 

(ii)  The  ASE  module  which  performs  the  new  features  overviewed  above  and  detailed 
below. 

(iii)  The  design  optimization  part  which  reads  the  ASE  constraints  and  their  sensitivities, 
combines  them  with  those  of  the  other  ASTROS  disciplines,  and  defines  design 
changes,  unless  the  process  is  converged.  This  part  will  remain  identical  to  the  current 
one,  except  for  some  possible  minor  changes. 

The  input,  output  and  methods  employed  in  the  blocks  of  the  ASE  module  and  the  new  bulk  data 

entries  for  each  block  are  described  below. 

4.4  Generalized  Matrices 

Generate  the  structural  and  aerodynamic  generalized  matrices  for  ASE  analysis. 

Input: 

(i)  g-set  stiffness  and  mass  sensitive  matrices  [DKVI]  and  {DMVI\,  associated  with  the 
stractural  global  design  variables. 

(ii)  The  current  values  of  the  global  design  variables  (v,). 

(iii)  Aerodynamic  AIC  matrices,  calculated  by  ZAERO  module  at  user  defined  reduced- 
frequency  values  (ki). 

(iv)  Kinematic  modes  generated  by  ZAERO  in  the  aerodynamic  grid: 

•  Control  modes  of  control  surface  unit  rotations. 

•  Gust  complex  modes  due  to  sinusoidal  gusts. 

•  Load  modes  for  calculating  section  loads  such  as  bending  moments  or  torques. 

•  Spline  matrices  between  the  aerodynamic  and/-set  stractural  coordinates. 

•  Stress/strain  coefficient  matrices. 


Output: 

The  coefficient  matrices  in  the  frequency-domain  open-loop  undamped  aeroelastic  equation  of 
motion  excited  by  control  surface  { 5c)  and  gust  inputs  [Wg } 

[  {icof  [Mhh]  +  [Khh]  -  q^lQhh  (i^)]  ]  (^  =  -  ( i^hc]  -  q^iQhc  0^)] )  (^d 

('lAi) 

^  0O 


and  the  aerodynamic  and  inertial  load  coefficient  matrices 
[Qui  (i^)] ,  [Qlc  0’^)] ,  [Qig  (f^)] ,  [Mit] ,  and 
This  block  expands  the  generalized  matrices  used  in  ASTROS  Version  11  (Ref  12). 
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Methods: 

The  ASE  module  will  allow  several  design  iterations  with  the  same  baseline  modes.  The 
operations  in  this  block  are  divided  into  two  cases: 


(i)  For  the  baseline  structure: 

•  Calculate  the  normal  modes  [(pafl  in  the  a-set  and  recover  to  the  /-set  [^]. 

•  Transform  [^^7,]  to  the  aerodynamic  grids  as  well  as  and  [(pki]  to  the  structural 
/-set  grids  [^J  and[0j. 

•  Use  AICs  and  k-set  modes  to  calculate  the  generalized  aerodynamic  force  matrices 
[2m]  >  [2ac]  >  \Qh^  >  [Qui]  >  [QiA  >  and  for  the  tabulated  k;  values. 

•  Recover  the  /-set  [^] ,  [^c]>  and  [^]  to  the  g-set. 

•  Calculate  sensor  modal  displacements  [^m]- 

•  Use  [DKVI\,  [DMVI\,  and  the  g-set  ,  [(pg^,  and  to  calculate  the  generalized 
stiffness  and  mass  sensitivity  matrices  associated  with  [Khi^,  [Mhh\,  [Mfic\,  [MlaJ*  and 

[Mi^],(e.g.[OTtt,.]=  M lDKVI][^i,]  ). 

•  Extract  the  diagonal  \Mhh\b  and  [ATaaIa  and  assemble  the  baseline  fMAclfe>  and 

[Mi^h  matrices,  where  subscript  b  represents  the  baseline  model. 


(ii)  For  a  modified  structure: 

Modify  the  generalized  stiffness  and  mass  matrices  according  to  the  generalized  sensitivity 
matrices  and  the  current  values  of  the  design  variables: 


e.g.  [KhhIm  =  [Khh\b  +  ^  (v.-  -  Vib )  [DKnh^  (4-4.2) 

where  the  subscript  M  denotes  the  modified  stracture. 

Presently  in  ASTROS,  {DKVI\  and  [DMVI\  are  not  constant  for  nonlinear  design  variables;  i.e. 
design  variables  that  cannot  be  factored  from  the  elemental  stiffness  and  mass  matrices  Kee  and 
Meg.  Under  the  new  paradigm,  we  can  still  make  use  of  the  ori^al  modes  to  reduce  the  current 
model,  but  we  must  augment  the  modal  stiffness  and  mass  updating  algorithm: 

[^gg]A/  ~  ^^eefijxd  -4  ^  P ij  Vj  +  AKee^  (4.4.3) 

j  »■ 

where  A  denotes  an  operator  for  matrix  assembly,  the  first  two  terms  represent  the  original 
(factored)  design  variables  stiffness  matrices  associated  with  the  nonlinear  design  variables  (Ref 
18),  and  the  Pij  represent  the  ASTROS  design  variable  linking  coefficients. 

During  the  first  design  iteration,  we  can  compute  and  save 


dKhh 


3vi 


\fact 


+  A  ES  Pij  ^^actl 


Then  the  [Khh]M  can  be  directly  computed  on  all  subsequent  design  iterations  as: 


(4.4.4) 
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[^Wi]m=  X 


(4.4.5) 


dKt^ 


dvi 


Vi  +  <j>gh\A  ]  ^gh 


jfact 


where  the  second  term  is  relatively  inexpensive  to  compute  since  AKeer,„d  typically  very 
sparse.  (j>gh  represents  the  original  modal  matrix  and  [A/aa]m  follows  in  an  identical  manner. 

To  complete  the  update  of  the  system,  new  first  order  modal  system  matrix  sensitivities  are 
need^.  In  the  absence  of  nonlinear  desi^  variables,  the  sensitivities  require  no  update  and  are 
identical  to  those  of  the  first  design  iteration.  However,  the  nonlinear  design  variables  introduce 
an  additional  term: 


dKhh' 

— 

dKhh 

T 

+  ^gh 

.  . 

M 

dvi 

fact 

L  ;  J 

<Pgh 


(4.4.6) 


whtTcdKee  /  is  the  finite-difference  sensitivity  of  the  elemental  stiffness  to  the  nonlinear,  local 

design  parameter,  U .  again  are  the  original  modal  matrix.  The  mass  sensitivities  follow 
directly. 

For  the  traditional  ASTROS  disciplines  STATICS  and  static  aeroelasticity  (SAERO),  the  current 
ASTROS  approaches  will  be  used.  The  processing  will  be  unaffected  by  the  modal  reduction 
scheme  used  in  the  ASE  analysis.  For  the  MODES  analysis,  the  reduced  Khh  and  ^hh  will  be 
used  to  compute  updated  eigenvalues  for  purposes  of  ASTROS  frequency  constraints.  If  no 
frequency  constraints  exist,  the  MODES  analyses  will  not  be  performed  after  the  first  design 
iteration. 

4.5  Control  Model 

Construct  the  s-domain  transfer  function  matrix  [T  (5)]  that  relates  actuator  outputs  to  sensor 
inputs: 

j5Wl=rr(j)l{yM)  f45n 


Input: 

(a)  Actuator  transfer  functions  of  the  form: 

_  _ ^io _ 

Uaci(s)  +  ai2  +  an  S  +  (4.5.2) 

The  order  difference  between  the  numerators  and  denominators  is  3  to  allow  for  the  appearance 

of  5n  Si  as  independent  states,  and  to  allow  for  direct  connection  to  acceleration  sensors. 

A  high  order  actuator  can  be  defined  by  adding  cascade  transfer  function  in  (c)  below. 

(b)  Sensor  transfer  function. 

(c)  Additional  filters  expressed  in  either  s-domain  transfer  function  or  state-space  form. 

Output: 

The  terms  of  [T  {s)],  each  expressed  as  a  series  of  s-domain  transfer  functions. 
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Method: 

Most  of  the  processing  is  just  relating  the  input  data  to  the  individual  terms  of  [T  {s)\-  State-space 
filter  components  are  converted  to  polynomial  form. 

New  Bulk  Data  Entries: 

The  introduction  of  the  control  model  requires  an  interface  between  the  ASE  module  and  the 
user.  Four  new  bulk  data  entries  are  defin^  and  described  as  follows: 

•  Input  data  entry:  SENSOIt 

Description:  Define  sensor  location  and  sensor  dynamics. 


Format  and  Exarrple 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

SENSOR 

IDSEN 

BID 

GD 

CD 

ITYPE 

lORD 

IDYNSET 

SENSOR 

1000 

1 

13 

3 

1 

0 

7 

IDSEN 

BID 

GD 

CD 

ITYPE 


lORD 


IDNYSET 


Sensor  ID 

Bomdary  condition  ID 
Grid  ID 
Component  ID 

Integer  indicating  whether  sensor  measures  linear  or  angular  motion 
=  0  determined  by  CD 
=  1  linear  motion 

-  2  angular  motion 

-  0  position  sensor 

-  1  rate  sensor 

=  2  acceleration  sensor 

^  0  DYNSET  set  ID  in  which  the  coefficients  of  the  numerator  and 
denominator  polynomial  are  specified 
=  0  perfect  sensor^  do  not  include  sensor  dynamics 


•  Input  data  entry:  DYNSET 

Description:  Define  the  values  of  the  numerator  and  the  denominator  polynomials  of 

transfer  function. 


Format  and  Example 


1 

3 

4 

5 

6 

7 

8 

9 

10 

DYNSET 

ID 

NUM 

NDE 

AN(1) 

AN(2} 

• . . 

.  • . 

. . . 

+ABC 

DYNSET 

31 

6 

8 

0.1 

2.0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DYNSET 

AN(NDM) 

AD(1) 

AD  (2) 

... 

... 

... 

AD  (NDE) 

+ABC 

1.9 

3.1 

0.0 

0.5 

ZisMi 

ID 

NUM 

NDE 

AN(i) 

AD(i) 


set  ID  of  DYNSET 
order  of  numerator 
order  of  denominator 

coefficient  of  polynomial  for  niomerator 
coefficient  of  polynomial  for  denominator 


Remarks: 

(1)  The  transfer  function  between  input  X  and  output  Y  is  defined  as: 


Y  AN(i)+AN{2)s  +—  +AN(NUM)s^VM-J  +sNUM 
^  ADil)+AD(2)s  +■■■  +AD{NDE)sNDE-1  +sNDE 


42 


STTjR  Final  Rjji05089S 


(2)  NUM  must  be  less  than  or  equal  to  NDE. 
•  Input  data  entry:  ACTOair 
Description:  Define  actuator  dynamics 

Format  and  Example 


ACTUAT 


ACTUAT 


Field: 

ID 

LABEL 

AD(1), 
AD(2), 
AD  (3) 


2 

3 

4 

5 

6 

ID 

LABEL 

AD(1) 

AD  (2) 

AD  (3) 

101 

ELEV 

0.0 

0.1 

0.5 

Actuator  ID 

Aerodynamic  control  surface  label  defined  in  AESURF  to  which  the  actuator  is 
attached 

Denominator  coefficients 


The  transfer  f motion  of  the  actuator  is  defined  as  _ 

^  AD{l)+AD{2)s  +AD{3)s2^s^ 

•  Input  data  entry:  FILTE31 

Description:  Define  type  of  filter  for  coirpensation  of  the  signals  from  the  sensors  prior 

to  sending  them  to  the  actuators. 

Format  and  Example 


2  3  4  5 


FILTER  I  IDFTR  |  ITYPE  A1  A2 


FILTER  I  10  I  3  I  10.0  1  0.01 


IDFTR  ID  of  FILTER 

iraPE  =  1  notch  filter:  - l+G^/Al^) - 

1+2s(A2/A1)+G2mi2) 

=  2  integral  filter:  Al/S 

—  3  lead-lag  (lag— lead)  filter:  ^ 

l+(A2)y 

<0  polynomial  filter  Y_ANil)+ANi2)s  -h-  +AN{NUM)sNUM-1  +  sNUM 

^  AD{1)  +AD  {2)s  +...+ AD  (A©£>M)£:-i  +  sNDE 

|lllfPE|  represents  the  set  ID  of  a  DYNSET  entry  which  defines  the  numerator  and 
denominator  polynomials  of  the  polynomial  filter 
Al,  A2  values  used  for  IDFTR  =1,  2,  or  3 

4.6  Rational  Aerodynamic  Approximation 

Approximates  the  combined  generalized  unsteady  aerodynamic  force  coefficient  matrix: 


by  a  rational  function  of  ik  which  has  the  form: 

[Q  (/^)]  =  [Ao]  +  [Aj]  (ik)  +  [A2]  {ikf  +  [D]  {[/]  {ik)  +  [E]  {ik) 


(4.6.1) 


(4.6.2) 
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where  the  coefficient  matrices  are  real. 


Input: 

(i)  Generalized  matrices  generated  in  block  2. 1 . 

(ii)  Actuator  transfer  function  from  block  2.2. 

(iii)  Type  of  data  weighting  and  design  flight  conditions  if  physical  weighting  is  used. 

(iv)  Approximation  order  and  constraints. 

(v)  The  diagonal  [/?] . 

(vi)  Initial  guess  of  [£>]  and  number  of  approximation  iterations. 


,  and  M  of  Eq  (4.6.2).  Partitions  of  [Aj]  are  as  in  Eq  (4.6.1), 


Output: 

(i)  [Ao],  [Ai],  [A2],  [E] 

\Dh] 

partitions  of  [O]  are  lAlJ  ,  and  partitions  of  [£]  are  [Eh  Ec  Eg] . 

(ii)  Table  of  maximum  weight  assigned  to  each  term  of  [Q  (i^)] . 


Method: 

The  minimum-state  rational  approximation  method  with  physical  weighting  and  constraints  is 
described  in  detail  in  (Refs  41  and  42).  The  approximation  process  is  based  on  a  series  of 
weighted  lease-square  solutions: 

(2  U  Ul)  =  1  [6l,  (4.6,3) 

where  (x* }  is  a  vector  of  unknown  coefficients,  [z?*]  is  based  on  the  aerodynamic  tabulated  data 
associated  with  ki,  and  [w*]/  is  a  diagonal  matrix  with  weights  associated  with  the  data  terms  in 
[b*]i.  When  no  constraints  are  applied,  the  process  starts  with  an  initial  [f>]  and  solves  [Ao],  [A;], 
[A2],  and[E],  column  by  column,  by  applying  Eq  (4.6.3)  with 

-k^I  kl[D\(k{\I\*[RfY 

{»'],=  I 

where  {Fj{ki)}  and  [Gj{Ki)}  are  real  and  imaginary  parts  of  the  j-th  column  of  the  tabulated 
matrix  [Q{iki)].  The  resulting  [£]  is  used  to  calculate  [Ao],  [A;],  [A2],  and  [f>],  row  by  row,  by 
solving  Eq  (4.6.2)  with 

I  0  -kfl  kf  E'^(kf{i\+[RfY 
0  kil  0  -ki  (kf  [/]+  [RfY  [R]  \  ’ 


[A*l  = 


{x*)  = 
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These  M  [£]  [£>]  iterations  are  repeated  until  convergence  is  obtained. 

The  users  will  have  the  option  to  use  Roger’s  approximation  which  can  be  posed  as  a  special 
case  of  Eq  (4.6.2),  see  (Ref  42). 

Each  column  of  [q  (i^/)]  will  be  approximated  with  up  to  3  approximation  constraints  selected 
from  the  following: 

(i)  Steady  aerodynamic  match  at  k  =  0. 

(ii)  Either  a  real-part  data-match  constraint  at  a  non-zero  it/  value  or  a  j^A2j]  =  0  constraint. 

(iii)  An  imaginary-part  data-match  constraint  at  a  nonzero  it/  value. 

There  will  be  three  weighting  options: 

(i)  Uniform  weighting,  *]/  =  [/]  (4.6.6) 

(ii)  Data-normalization  weighting  Wy.  = - ^  ,  ,, — r  (4.6.7) 

maxi{\Qij{iki)\,e} 

where  e  is  a  user-defined  small  positive  parameter. 

(iii)  Physical  weighting,  as  detailed  in  (Refs  41, 47). 

When  option  (iii)  is  selected,  the  weights  of  the  h  column  partition  in  Eq  (4.6.1)  are  based  on 
derivatives  of  the  system  matrix  in  the  left  hand  side  of  Eq  (4.4.1).  The  weights  of  the  c  column 
partition  of  Eq  (4.6.1)  are  based  on  Nyquist  return  signals  with  the  control  system  assumed  to  be 
represented  by  the  actuator  transfer  functions,  Eq  (4.5.2).  The  weights  assigned  to  the  gust 
coluxnns  are  based  on  a  selected  response  to  Dryden’s  gust  spectrum. 

The  physical  weighting  assigns  measures  of  aeroelastic  importance  that  can  be  used  to  select 
modes  for  residualization  or  trancation. 

New  Bulk  Data  Entries: 

One  new  bulk  data  entry  is  required  to  define  the  input  parameters  for  the  minimum  state 
technique  and  is  described  as  follows: 

•  Input  data  ontzy:  KIST 

Description:  Define  parameters  for  rational  approximation  of  unsteady  aerodynamic  forces 

by  minimum  state  technique. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

MIST 

ID 

KLIST 

M 

NKF 

ITMAX 

IWE 

WCUT 

NWD 

+  ABC 

MIST 

1 

2 

3 

1 
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Eioldi. 

ID  ID  nijirriber 

KLIST  ID  number  of  AEFACT  entry  specifying  a  list  of  hard  point  reduced  frequencies 

M  Number  of  approximation  roots 

NKF  =  1  all  real  and  imaginary  parts  of  generalized  forces  are  matched  at  the 

last  hard  point  reduced  frequency  list  in  FEFACT  entry 
>  1  all  real  and  imaginary  parts  matched  at  the  NKF'th  reduced 
frequency  listed  in  AEFACT  entry 

=  0  real  parts  are  matched  at  the  last  hard  point  reduced  frequency  and 
imaginary  parts  at  K(2) 

ITMAX  <  0  [^21=0,  imaginary  parts  are  matched  at  k  (-NKF)  number  of  D  E  — >  D 

iterations 

IWE  =  1  use  weight  matrices  from  a  previous  run 

=  2  calculate  weight  matrices  based  on  physical  weighting 
=  3  calculate  weight  matrices  by  data-normalization 
WCUT  Minimal  maximum  absolute  value  of  weighted  aerodynamic  term  (0.0  <  WCUT  <  1.0) 

NWD  Number  of  weight --peak  widening  cycle 

R(i)  M  number  of  approximated  roots  (distinct  negative  values) 


Calculates  gain  margins,  phase  margins,  and  singular  values  of  the  MIMO  aeroservoelastic 
system. 


Input: 

(i)  Retrieve  generalized  structural  matrices  of  Eq  (4,4.1)  from  block  2.1. 

(ii)  Control  transfer  function  from  block  2.2. 

(iii)  Modal  displacements  of  sensor  inputs 

(iv)  Definition  of  the  control  design  variables. 

(v)  Rational  approximation  of  [<2  (t^)] ,  Eq  (4,6.2)  from  block  2.3. 


Output: 

(i)  Gain  margins,  phase  margins  and  minimum  singular  values  in  the  user-defined  frequency 
range. 

(ii)  Derivatives  of  these  control  margins  with  respect  to  structural  and  control  design  variables. 


Method: 

The  aeroelastic  plant  in  Eq  (4.4.1)  is  expanded  to  include  the  diagonal  matrix  of  actuator  transfer 
functions  [TacC^)]  and  the  diagonal  matrix  of  sensor  transfer  functions  [T,(5)].  The  resulting 
transfer  function  of  the  plant,  from  actuator  input  to  sensor  output  is 

[G  (s)]  =  [T.  (sj  [fa]  [Cm  (sr  («_[efc(s)  -  [Mte]]  [r^  (sj)  (4.7.1) 

The  remaining  control  feedback  system  \k  (s)],  which  relates  actuator  inputs  { Mac}  to  command 
error  {e}  is  based  on  the  filters  defined  in  block  2.2.  The  MIMO  control  aeroservoelastic  system 
block  diagram  is  given  in  Figure  4.7.1. 
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Figure  4.7.1  Block  Diagram  of  the  MIMO  ASE  System. 

The  command  errors  are  defined  as  {e}  =  {}’coot}  -  {y}-  To  analyze  the  system  stability,  it  is 
assumed  that  {ycom}  =  0  such  that  the  input  return  signal  is  {uac}  =  -  M  [G]  {u}  and  the  input 
return-difference  matrix  is  [/  +  KG],  As  suggested  by  Ref  48,  the  minimum  singular  values  of 
[I  +  KG],  denoted  by  <t  {I  +  KG),  and  those  of  the  inverse  return-difference  matrix 

O’  (/+  {KGy^),  will  be  used  as  measures  of  guaranteed  simultaneous  stability  margin.  The  s 
values  will  be  calculated  over  a  range  of  s  =  ict)  values.  Singular-value  constraints  at  several 
frequency  points  and  calculated  derivatives  of  o  (ito)  at  these  frequencies  with  respect  to  the 
stmctural  and  control  design  variables  will  be  included  in  the  optimization  process. 

To  conform  with  more  classical  control  design  requirements,  gain  and  phase  margins  and  their 
sensitivities  will  be  calculated  by  Nyquist-type  andyses  with  die  control  loops  opened  one  at  a 
time  at  the  actuator  inputs. 

New  Bulk  Data  Entries: 

Block  2.4  allows  the  user  to  define  constraints  of  the  ASE  system  in  terms  of  gain  margins,  phase 
margins,  and  singular  values.  Three  new  bulk  data  entries  are  defined  and  described  as  follows: 

•  Input  data  entzy:  DOONMAR 

Description:  Define  the  gain  margin  constraint  of  the  ASE  system  with  the  control  loops 

opened  one  at  a  time  at  the  specified  actuator  inputs. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DCONMAR 

CTSET 

crypE 

IDACT 

RMAR 

IFREQ 

DCONMAR 

1 

UPPER 

31 

3.0 

10 

Eield:. 

CrSET 

CTYPE 

IDACT 

KMAR 

IFREQ 


Constraint  set  identification 

Constant  type^  either  UPPER  for  upper  bound  or  LOWER  for  lower  bound  (Text;, 
default  -  LOWER) 

Actuator  ID 

Required  gain  margin  value  (dB) 

Frequency  set  identification  number  of  bulk  data  FREQ,  FREQl,  or  FREQ2 


Remarks: 

The  frequencies  listed  in  IFREQ  represent  the  frequency  points  at  which  the  gain  margins 
cannot  exceed  the  constraint  value  RMAR. 


•  Input  data  entry:  DCONPHS 

Description:  Define  the  phase  margin  constraint  of  the  ASE  system  with  the  control  loops 

opened  one  at  a  time  at  the  specified  actuator  inputs. 

Format  and  Exanple 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DCONPHS 

CTSET 

CTYPE 

IDACT 

RPHS 

IFREQ 

DCONPHS 

2 

UPPER 

31 

30.0 

10 

El&lsh. 

CTSET 

CTYPE 

IDACT 

RPHS 

IFREQ 


Constraint  set  identification 

Constraint  type  either  UPPER  for  upper  bound  or  LOVJER  for  lower  bound.  (Text, 
Default  =  Lower) 

Actuator  ID 

Required  phase  margin  value  (degrees) 

Frequency  set  identification  number  of  bulk  data  FREQ,  FREQl,  or  FREQ2 


Remarks:, 

The  frequencies  listed  in  IFREQ  represent  the  frequency  points  at  which  the  phase  margins 
cannot  exceed  the  constraint  value  RPHS. 


•  Ixxput  data  entry:  DOQNSIN 

Description:  Define  the  singular  value  constraint  of  the  ASE  system  in  a  specified 

frequency  range. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DCONSIN 

CTSET 

CTYPE 

SING 

IFREQ 

DCONSIN 

1 

0 

0.3 

10 

Field: 

CTSET  Constraint  set  identification 

CTYPE  =  0  singular  value  of  return-difference  matrix 

=  1  singular  value  of  inverse  return-difference  matrix 
SING  Required  minimum  singular  value  (positive) 

IFREQ  Frequency  set  identification  number  of  bulk  data  entry  FREQ,  FREQl, 

or  FREQ2 

The  frequency  listed  in  IFREQ  represent  the  frequency  points  at  which  the  frequency  points  at 
which  the  singular  value  must  be  greater  than  SING. 

4.8  state  -  Space  ASE  Model 

Constructs  the  ASE  state-space  model  for  flutter  and  gust-response  analyses  and  export  for 
control  analyses. 

Input: 

(i)  Generalized  structural  matrices  of  Eq  (4.4. 1)  from  block  2. 1 . 

(ii)  Control  transfer  functions  and/or  state-space  matrices  from  block  2.2. 

(iii)  Rational  approximation  of  \Q-  (^^H,  Eq  (4.6.2),  from  block  2.3. 

Output: 

(i)  State-space  l^ae],  [Bae],  and  matrices  describing  the  aeroelastic  system  (plant  + 
actuators  +  sensors). 

(ii)  State-space  \Ac] ,  [^c] ,  [Q] ,  and  matrices  describing  the  control  system. 

(iii)  Closed-loop  control-augmented  system  matrix  H . 

(iv)  Output  matrix  for  gust  response  analysis. 

(v)  Sensitivities  of  the  above  matrices. 

Method: 
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The  open-lcx)p  state-space  equation  of  motion  of  the  aeroelastic  system  excited  by  a  control 
surface  is 


[xp  I  -  [Ap]  {xp }  +  [Bp]  (Mp) 


where 


{^p)=  t 


r 

W=  6 


(4.8.1) 


[Ap]  -  -  +  qAhho \  -  Bhh  +  ,  -qoo\M\'^ 

^  OP 


[Bp\=  -Qo^lWlAhco]  .  -"^[WAhcll  ,  -[W  Mhc+^^[Ahc2] 

L  J 


,  M=M  +  ^[Am2]  „  /A..  ,  ^  ^  , 

•  Here,  \s/  is  the  vector  of  generalized  structural  displacements 
and  (b)  is  the  vector  of  control  surface  command  deflections.  The  plant  output  equation,  based 
on  the  sensor  modal  deflections  [^5/,]  is  described  in  the  form: 

{yp }  =  [Cp]  W }  +  [Dp]  [up }  (4.8.2) 

where  [Cp]  and  [Dp]  depend  on  the  type  of  the  sensor  (displacement,  velocity,  or  acceleration). 
The  actuator  of  Eq  (4.5.2)  is  modeled  in  state-space  form  as 


0 

1 

0  ' 

0 

0 

1 

0  ' 

-  “^io 

‘dil 

-a-a  . 

Wio] 

(4.8.3) 


where 
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Several  actuators  can  be  grouped  together  for  the  general  form 


—  [^ac]  {^ac}  [®acl  {^ac} 

Augmentation  to  the  plant  Eqs  (4.8.1)  and  (4.8.2)  yields 


xp  L 

\Ap\  [Bp] 

1  !+ 

■  0  ■ 

Xac  j 

0  [A^]J 

\  Xac  i 

.[Bac]. 

since  [xac]  =  {up},  and  the  output  (sensor  measurement)  equation  becomes 

(^l=[[Cp][Dp]]{  I  j 

A  general  transfer  function  of  a  control  component  is 

(^)  ^bos'^  +  +  •  ■  •+  bn 

^  Uc{s)  +  +  •••+an 

The  controller  canonical  form  realization  of  this  transfer  function  is 
{jCc}  =  M{jCc}  +  {Bc}  {Uc} 
yc  =  [Cc](:tc}  +  [Dc]{uc} 

where 


■  0  1 

0  ■ 

("1 

[Ac]  = 

0  '  6 
.  -  On  -  dn-l  •  •  • 

1 

1 

:  \ 

[Cc]  =  [{t>n  -  boUn)  [bn-i  -  boUn-l  )•••(*/-  boai)[ 

Dc  =  bo 

Application  of  Eq  (4.8.8)  to  the  sensor  transfer  functions  and  augmentation  with  Eq 
the  state  space  equations  of  the  aeroelastic  system 

{jCfle}  =  [Afle]  "I"  [BqJ  (Wae) 

{yae)  —  [CflJ  {Xae} 

Finally,  the  state-space  equations  of  the  control  filters  can  be  grouped  in  the  form 


(4.8.4) 

(4.8.5) 

(4.8.6) 

(4.8.7) 

(4.8.8) 

.8.5)  yield 

(4.8.9) 
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{jCc}  ~  D^c]  {-^c}  [^c]  {^c} 

{}'c)=[Cc]{Xc)+[£>c]{Mc} 


(4.8.10) 

Augmentation  of  Eq  (4.8.10)  to  Eq  (4.8.9)  with  {uc}  =  {yae}  yields  the  general  open-loop  control 
form 

{i}  =  [A]  +  {«} 

(4.8.11) 

Connection  of  the  inputs  (m)  to  the  outputs  {y }  through  a  gain  matrix,  {u}  =  [G]  (y  ),  yields  the 
closed-loop  equation 

{i}  =  [A]{x}  (4  3  12) 


where 


H=W+[B][G][q 

The  block  2.5  assembly  of  the  state-space  ASE  system  matrix  sensitivities  with  respect  to 
structural  and  controller  design  parameters  must  be  evaluated  if  any  constraints  associated  with 
the  ASE  analyses  are  retained  as  active.  With  the  advent  of  design  variables  associated  with  the 
control  law,  the  derivatives  must  account  for  both  FE  matrix  perturbations  (which  are  well 
understood  and  have  been  extended  to  support  the  iterative  formation  of  a  generalized  modal  and 
its  first  order  derivative)  and  control  parameter  perturbations.  The  control  parameters  include 
gains  and  the  coefficients  of  the  polynomials  for  each  acmator/sensor.  The  user  will  describe 
which  of  these  parameters  is  available  for  design  and  its  range  of  values  [v^,  v,^. 

The  sensitivities  of  the  ASE  constraints  (or  the  objective  function)  will  be  solved  analytically 
based  on  the  sensitivities  of  the  state-space  [A^J,  and  matrices  describing  the  system 
and  of  the  state-space  controller  matrices  [AJ,  [5J,  [CJ  and  [DJ.  It  will  be  assumed  that  the 
structural  design  parameters  affect  only  [1^  and  [M]  and  the  controller  design  parameters  do 
NOT  affect  [i^]  and  [M]  or  the  rational  aerodynamic  approximation.  Traditional  ASTROS 
constraints  wM  not  be  affected  except  that  the  vector  of  design  variables  must  be  kept  separated 
in  the  structural  and  controller  parameters. 


The  cross-sensitivities  will  be  ignored  for  purposes  of  sensitivity  calcidation. 
New  Bulk  Data  Entries: 


•  Input  data  entry:  GZUU 

Description:  Define  gain  of  a  control-sensor  pair. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

GAIN 

ID 

LABEL 

IDSEN 

GN 

PHASE 

IDFTR 

GAIN 

4 

ELEV 

20 

0.7 

0 

10 
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Eig.ldi 

CrSET  ID  of  gain 

LABEL  Alphanumeric  string  for  aerodynamic  control  surface  defined  in  input  data  entry 

AESURF 

IDSEN  Sensor  ID  number  defined  in  input  data  entry  SENSOR 

GN  Value  of  gain 

PHASE  Phase  error  (degrees)  in  transfer  fimction 

IDFTR  ID  of  filter  associated  with  the  current  control-sensor  pair 


•  Input  data  entry:  DEISGAIN 

Description:  Define  gain  value  of  a  control-sensor  pair  as  a  design  variable. 


Format  and  Exairple 


1 

3 

4 

5 

6 

7 

8 

9 

10 

DESGAIN 

DVID 

GID 

VMIN 

VMAX 

VINIT 

LABEL 

DESGAIN 

3 

4 

0.1 

0.3 

0.2 

GAINl 

FifflA; 

DVID  Design  variable  identification 

GID  ID  of  data  entry  GAIN 

VMIN  Minimum  allowable  value  of  the  design  variable 

VMAX  Maximum  allowable  value  of  the  design  variable 

VINIT  Initial  value  of  the  design  variable  (VMIN  <  VINIT  <  VMAX) 

LABEL  Optional  user-supplied  label  to  define  the  design  variable  (text) 

•  Input  data  entry:  DESOOKF 

Description:  Define  a  coefficient  of  the  polynomials  associated  with  the  transfer  function 

of  a  control  filter  as  a  design  variable. 


Format  and  Exaitple 


1 

3 

4 

5 

6 

7 

8 

9 

10 

DESCOEF 

DVID 

IDFTR 

ITYPE 

N 

VMCN 

VMAX 

VINIT 

LABEL 

DESCOEF 

1 

2 

1 

3 

0.1 

0.3 

0.2 

COEl 

Field: 

DVID  Design  variable  identification 

IDFTR  ID  of  bulk  data  entry  FILTER 

ITYPE  =  1  define  the  coefficient  A1  as  design  variable  for  notch,  integral, 

proportional  plus  derivative,  or  lead-lag  filter 
=  2  same  as  ITYPE  =  1  but  for  A2 

-  3  define  one  of  the  coefficients  of  the  numerator  polynomial  as  a 
design  variable 

“  4  same  as  ITYPE  =  3  but  for  the  denominator  polynomial 
N  The  N-th  coefficient  of  the  numerator  or  denominator  (Used  only  for 

ITYPE  =  3  or  4) 

VMIN  Minimum  allowable  value  of  the  design  variable 

VMAX  Maximum  allowable  value  of  the  design  variable 

VTNIT  Initial  value  of  the  design  variable  (VMIN<VINIT<VMAX) 

LABEL  Optional  user-supplied  label  to  define  the  design  variable  (text) 


4.9  Flutter 

Calculates  flutter  velocity  and  the  same  flutter  constraints  as  in  ASTROS  Version  11,  but  based 
on  a  state-space  formulation,  with  or  without  a  control  system. 

Input: 

Closed-loop  system  matrix  from  block  2.5. 

Output: 
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Flutter  speed. 

Damping  constraints  and  their  sensitivities. 

Method: 

The  flutter  analysis  is  based  on  the  eigenvalues  of  the  closed-loop  system  matrix  of  Eq 
(4.8.12)  or  the  open-loop  matrix  [Ap]  of  Eq  (4.8.1)  when  there  is  no  control  system.  Root  locus 
analysis  is  performed  for  the  flutter  velocity  with  fixed  density  or  for  the  flutter  density  with 
fixed  velocity.  The  definition  of  the  damping  constraints  and  the  sensitivities  of  the  active 

constraints  will  be  similar  to  those  in  Version  1 1,  based  on  the  eigenvectors  of  or  [Ap]  at  the 
required  constraint  points. 

4.10  Gust  Model 

Constructs  a  state-space  model  where  the  input  signal  represents  a  white  noise  process  and  the 
output  signal  represents  Dryden’s  or  von  Karman's  continuous  gust  velocity. 

Input: 

-  User  input  parameters  for  either  Dryden’s  or  von  Karman's  power  spectral  density  (PSD) 
functions. 


Output: 

State-space  gust  matrices  [Ag],  \Bg\,  and  [Cg]. 

Method: 

The  construction  of  the  gust  model  is  based  on  (Ref  49).  Either  Dryden's  or  von  Karman's  gust 
PSD  functions  are  represented  by  a  Laplace  domain  transfer  function,  from  white  noise  signal  w 
to  gust  velocity  Wg.  The  order  of  the  transfer  function  denominator  is  larger  than  that  of  the 
numerator.  The  resulting  gust  model  has  the  form 


{%)  “  {^}  w 

with  three  states  for  Dryden's  gust  and  four  states  for  von  Karman’s  gust. 
New  Bulk  Data  Entries: 

One  new  bulk  data  entry  is  defined  for  the  continuous  gust  input. 


(4.10.1) 


•  Input  data  entry:  GDSTl 

Description:  Defines  Dryden's  or  von  Karman's  gust  PSD  functions  for  ASE  gust  analysis. 


Format  and  Exairple 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

GUSTl 

SID 

SYMXY 

ITYPE 

MIST 

WG 

XO 

V 

GUSTl 

1 

1 

1 

10 

1.0 

10.0 

1000.0 

SID  Gust  set  identification  number 

SYMXY  =  1  vertical  gust 

-  2  lateral  gust 

ITYPE  =  1  von  Karmon's  gust  PSD  function 

=  2  Dryden's  gust  PSD  function 

MIST  ID  of  bulk  data  entry  MIST  for  rational  approximator  of  the  harmonic 

gust  forces 


53 


STTR  Final  Rpt050896 


WG 

xo 

V 


Scale  factor  of  the  gust  velocity 

Location  of  reference  plane  in  aerodynamic  coordinates 
Velocity  of  vehicle 


4.11  Gust  Response 

Calculates  the  Root-Mean- Square  (RMS)  values  of  gust  response  parameters  and  their 
sensitivities. 

Input: 

-  Closed-loop  system  matrix  from  block  2.5. 

-  State-space  gust  matrices  from  block  2.7. 

-  Ration^  approximation  matrices  from  block  2.3. 

-  Sensitivities  of  the  system  matrices  from  block  2.5. 


Output: 

-  Gust  response  constraints  and  their  sensitivities. 

Method: 

The  gust  matrices  of  Eq  (4.10.1)  are  augmented  to  either  the  plant  matrix  [A^]  of  Eq  (4.8.1)  or  the 
closed-loop  matrix  [A]  of  Eq  (4.8.12),  as  detailed  in  (Ref  8).  The  resulting  system  has  the  form 

U)=WI^I+(-8w!»>' 

tyl=HW 

The  root-mean  square  values  of  the  response  parameters  are  based  on  the  state  covariance  matrix 
[X]  =^{x}  (x}^].  When  Eq  (4.11.1)  is  excited  by  a  unit  intensity  white  noise,  [X]  satisfies  the 
matrix  Lyapunov  equation 


[A][X]  +  [X][Ar=-{Bw)(fiwf 

The  covariance  matrix  associated  with  {y}  of  Eq  (4.11.1)  is 

m^£[{)>)(yF]=[Cl[X][Cf 


(4.11.2) 


(4.11.3) 


The  diagonal  of  [T]  contains  the  mean-square  output  response  parameters  of.  The  differentiation 

dlB  ) 

of  Eq  (4.11.2)  with  respect  to  a  design  variable  Vj,  using  ■  ■  ---  =  0,  yields  another  Lyapunov 

dvi 

equation 

[A]  ^  +  ^[Af  =  --^[X]  -  [X]  (4.11.4) 

dvi  dvi  dvi  dvi 


which  can  be  solved  for 


d[X]^ 

9vi 


A  new  solution  is  required  for  each  design  variable,  but  they  can 


all  be  based  on  the  same  decomposition  used  in  solving  for  [X]  by  Eq  (4.11.2).  The 
differentiation  of  Eq  (4.11.3)  yields  the  response  sensitivity  derivatives 
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(4.11.5) 


^=[C\  ^[cY +^[x]  [CY + [X]  [X] 

dvi  dvi  dvi  dvi 


where  the  diagonal  terms 


are  usually  the  only  derivatives  of  interest. 


New  Bulk  Data  Entries: 

Two  new  bulk  data  entries  are  introduced  for  the  gust  response  constraints  of  the  RMS  value  of  a 
group  of  aerodynamic  components  and  the  RMS  displacements,  velocities  and  accelerations  of 
grid  point. 

•  Ixxput  data  entry:  DOONACD 

Description:  Defines  the  gust  response  constraint  of  a  group  of  aerodynamic  components  in 

terms  of  forces  and  moments. 


Format  and  Exanple 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DCONACD 

SID 

lACORD 

ITYPE 

RMS 

IGUST 

DCONACD 

3 

10 

1 

100. 

1 

Eialdx 

SID 

lACORD 

ITYPE 

RMS 
IGDST 

Remarks; 

The  moment  axis  is  defined  by  XMCNT,-  YMCOT,  and  ZMCNT  in  ACORD  data  entry. 

•  Input  data  entry:  IXXXJBMS 

Description:  Defines  the  displacement,  velocity  or  acceleration  constraint  in  terms  of  RMS 

value  in  response  to  a  gust  input. 


Format  and  Example 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

DCONRMS 

SID 

DTYPE 

IGUST 

GID 

C 

DALL 

CTYPE 

DCONRMS 

1 

1 

1 

73 

1 

2.0 

UPPER 

Constraint  set  identification  number 

Identification  number  of  the  ACORD  bulk  data  entry  in  which  a  group  of 
aerodynamic  corrponents  is  defined 
=  1,  2,  or  3  :  force  along  x,  y,  or  z  direction 
-  4,  5,  or  6  :  moment  about  x,  y,  or  z  axis 
Constraint  in  term  of  RMS  value 

Identification  nxiniber  of  GUSTl  bulk  data  entry  where  the  gust  input  is  defined 


Field: 

SID  Constraint  set  identification  numbers 

DTYPE  —  1  displacement  constraint 

=  2  velocity  constraint 
=  3  acceleration  constraint 

IGUST  Identification  nunber  of  GUSTl  bulk  data  entry  where  the  gust  input  is  defined 

GID  Grid  ID 

C  Coirponent  ID 

DAIiL  Allowable  value 

CTYPE  Either  UPPER  or  LOWER  bound  (Text,  default  =  UPPER) 
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SECTION  5 


MODIFICATION  OF  ASTROS  FOR  ZAERO  AND  ASE  INTEGRATION 
5.1  Integration  Overview 

The  purpose  of  the  ASE  capabilities  in  ASTROS  is  to  facilitate  the  stractural  optimization  of  an 
aeroservoelastic  system  with  respect  to  objective  and  constraint  functions  associated  with  ASE 
Stability  and  Response  (Gust)  analyses.  The  parameters  used  to  effect  improvements  will  be  the 
usual  ASTROS  structural  parameters  and  with  the  addition  of  gains  and  polynomial  coefficients 
of  the  control  system.  The  control  system  architecture  is  fixed. 

It  is  important  to  note,  however,  that  the  proposed  ASE  capabilities  extend  beyond  the  straight 
forward  (black-box)  implementation  of  automated  design  subject  to  constraints  on  ASE  analysis 
responses.  Rather,  the  proposed  enhancements  are  intended  to  make  ASTROS  useful  in  the 
process  of  designing  a  ffight  vehicle  in  which  the  designer  wants  to  find  the  “best”  design  by 
exploiting  the  complex  synergies  between  the  structural  behavior  and  the  MIMO  control  system. 

Out  of  necessity,  this  is  an  interactive  process  that  is  dominated  by  perturbation  of  mathematical 
models  describe  by  parameters  that  cannot  be  updated  automatically.  For  example,  the  sensor 
and  acmator  locations,  filters  and  other  components  of  the  control  laws  do  not  readily  lend 
themselves  to  the  automated  design  techniques  of  traditional  structural  optimization.  The 
proposed  extension  of  ASTROS  wifi  enable  the  design  team  to  assemble  ASE  models,  optimize 
them  with  respect  to  stability,  control,  and  strength  requirements  and  to  compare  the  optimal 
designs  to  one  another  to  gain  insight  into  the  system’s  characteristics.  Because  many 
components  of  the  ASE  models  can  be  shared  and  are  computationally  intensive  to  produce, 
many  of  the  proposed  enhancements  deal  with  their  creation,  assembly,  management,  and  use  of 
these  components. 

The  proposed  ASE  analysis  approach,  minimum  state  rational  approximation  based  on  modal 
system  matrices,  requires  an  initialization  phase  in  ASTROS  to  create  the  modal  system 
matrices.  The  rational  approximation  and  the  generation  of  the  state-space  model  of  the  ASE 
system  will  then  be  performed  for  each  analysis  case.  The  controller  and  the  aerodynamics  may 
be  diffCTent  from  case  to  case;  e.g.  symmetric  flutter  at  M  =  0.9  in  one  case  and  antisymmetric 
gust  response  at  M  =  1.5  in  another.  Each  of  the  analysis  cases  could  have  different  control  laws 
as  well:  to  handle,  for  example,  separate  low-^  and  high-^  systems. 

Since  the  generation  of  the  ASE  model  for  a  given  analysis  case  is  computationally  intensive  and 
may  call  for  user  intervention,  a  requirement  of  this  integration  activity  is  that  ASTROS  produce 
“archival”  models  at  each  stage  of  the  process  for 

1)  subsequent  re-use  in  ASTROS  ASE 

2)  export  to  other  state-space  ASE  tools,  e.g..,  MATLAB  and  MATRIX  x. 

As  part  of  the  improvements  to  generalized  unsteady  aerodynamic  influence  coefficients  in  the 
ZAERO  module,  ASTROS  will  support  the  generation  and  archiving  of  a  very  general  set  of 
aerodynamics  that  will  provide  the  ASE  analyses  with  their  required  aerodynamic  influence 
coefficients.  The  other  remaining  activities  to  support  the  assembly  of  an  ASE  module  are  the 
creation  and  archival  of  modal  representations  of  the  basic  structurfi  system  matrices,  Khh,  Mhh, 
Bhk  and  of  the  response  matrices,  s  =  S/, «/,;  (where  5*  is  the  modal  stress  or  strain  displacement) 
and  the  generation  of  a  rational  aerodynamic  approximation  in  the  s-domain.  Each  of  these  two 
steps  have  implications  with  respect  to  the  current  ASTROS  paradigm. 
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To  create  the  rational  aerodynamic  approximation,  the  raw  AIC  data  for  a  given  set  {Mach 
number,  symmetry,  method,  modal}  need  to  be  fit  with  the  help  of  user  input  to  control  the 
fitting  process.  The  input  to  this  process  includes: 

1)  Mu,„  Kui,  and  Bu, 

2)  Qmi  (ik),  aerodynamic  geometry  &  control  surfaces 

3)  control  system  definition 

4)  algorithm  control  parameters 

As  part  of  this  process,  the  structural  equations  of  motion  for  the  open-loop  system  must  be 
assembled  and  solved.  Also,  at  this  time,  the  archived  AIC  matrix,  die  archived  aerodynamic 
geometry  and  control  surfaces  along  with  the  particular  disciplinary  parameters  describing  the 

analysis  subcase  (flight  condition,  gust  form  and  reference  data,  etc.)  are  used  to  produce  the  <{)kc, 

(t>Lg  (i^)  and  as  needed.  This  process  is  not  computationally  intensive  and  is  a  function  of 
parameters  which  the  user  may  want  to  vary  as  part  of  a  trade  study.  Thus,  these  data  are  NOT 
archived  as  part  of  the  aerodynamic  model.  The  results  are  a  set  of  coefficient  matrices  fitting 
the  aerodynamics. 

The  quality  of  this  fit  may  need  to  be  checked  for  the  set  of  responses  that  are  deemed  important 
in  the  current  analysis.  This  activity  may  need  to  be  performed  off-line  and  the  fit  may  need  to 
be  modified.  Finally,  the  fitting  coefficient  matrices  may  be  archived  for  reuse.  They  constitute 
a  particular  system  representation  that  includes 

1)  structural  geometry  and  element  properties  &  design  variables  “the  FE  model” 

2)  AIC  matrix  for  a  particular  aerodynamic  model  at  a  particular  (Mach,  number,  symmetry, 
method,  modal}. 

3)  a  spline  model  that  connects  1)  and  2) 

4)  a  particular  control  law  that  includes 

-  actuators 

-  sensors 

-  filters 

This  archived  approximation,  while  computed  based  on  the  current  design  point,  will  be  re-used 
during  the  subsequent  ASTROS  design  iterations,  with  the  user  free  to  update  the  fit  at  any  point. 

The  final  activity  to  complete  the  assembly  of  the  ASE  module  is  the  generation  and  archiving  of 
the  modal  representations  of  the  FE  model.  The  typical  modal  parameters:  [x], 

[Khiil,  and  [Bhhl  are  straight  forward.  However,  they  introduce  some  complexity  in  that  the 

particular  modes  [^g/,]  may  need  to  be  extracted  from  a  larger  archived  pool  of  modes  that  have 
been  computed.  They  do  not  necessarily  represent  all  the  modes  in  the  frequency  range 

[Xmin,  Xmal  from  the  archived  modal  analysis.  In  other  words,  to  support  ASE,  ASTROS  will 
allow  the  archival  of 


[Kiil  [Mii]  and[B,7] 


and  the  associated  [Ai,  A2---A,]  that  are  produced  by  normal  modes  analysis.  They  are 
associated  with  a  particular  structural  FE  model  and  a  particular  boundary  condition  (symmetric, 
antisymmetric).  The  ASE  disciplines  will  be  developed  such  that  they  “import”  these  data  for 
use.  The  assembly  of  the  h-sct  matrices  for  ASE  is  then  the  process  of  extracting  the  desirable 
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subset  of  [a,  •  ■  -Aj]  and  the  augmentation  of  the  model  with  additional  generalized  coordinates 
(extra  points,  residual  vectors,  etc.). 

5.2  ASE  Design  Overview 

In  design  optimization,  the  concept  of  re-using  the  archived  modal  representation  of  the  FE 
model  tecomes  more  complex.  The  aerodynamic  model  is  completely  reusable,  apart  from  some 
considerations  with  respect  to  the  rational  approximation  and  the  generalized  aerodynamic 
model.  Also,  the  generation  and  use  of  the  approximate  optimization  problem  in  ASTROS  has 
some  implication  with  respect  to  the  implementation  of  Blocks  2.1, 2.5, 2.6  and  2.8. 

5.3  Solution  Control 

ASTROS'  solution  control  will  be  enhanced  in  a  number  of  ways  to  support  the  proposed  ASE 
features.  These  enhancements  fall  into  seven  categories  (in  order  of  dependence): 


1)  Definition  of  one  or  more  unsteady  aerodynamic  cases 

2)  Archival/Import  of  aerodynamic  models 

3)  Definition  of  one  or  more  rational  aerodynamic  approximations  (RAA) 

4)  Archival/Import  of  RAA  data 

5)  Definition  of  one  or  more  control  margin  analyses 

6)  Definition  of  one  or  more  state-space  flutter  analyses 

7)  Definition  of  one  or  more  state-space  gust  analyses 


There  will  be  five  new  analysis  disciplines  to  create  aerodynamic  data  or  to  use  aerodynamics  to 
perform  response  analyses: 


ZAERO: 


MIST: 


MARGIN 


ASEFLUT: 


GUST: 


generates  \Ajj\  matrices  and  geometry  for  a  given  method  and  model  and  for  a 
series  of  {M,  k,  symmetry}  triplets  where  the  model  is  a  named  collection  of 
CAEROibulk  data  representing  a  configuration  and  the  method  is  selected  from 
I  ZONA6  CPM 

)  ZONA7  DLM 

j  ZTAIC 
\ZONA51U 

Minimum  state  approximation  to  fit  aerodynamic  data  for  one  or  more  Mach 
numbers  of  a  ZAERO  case. 

-  This  analysis  discipline  either  imports  an  archived  ZAERO  case  or  refers  to  a 
previously  defined  case  in  the  current  run. 

-  Imports  or  refers  to  a  MODES  case  for  structural  modes. 

Computes  gain  and  phase  margins  and  singular  values  of  a  MIMO  system,  for  a 
given  set  of  Mach  number,  boundary  conditions,  method,  and  model. 

-  This  analysis  discipline  either  imports  an  archived  MIST  case  or  refers  to  a 
previously  defined  MIST  case  in  the  current  form. 

-  In  optimization,  constraints  (upper  and  lower  bounds)  may  be  placed  on  some 
or  all  the  margins  and  singular  values. 

Performs  state  space  flutter  analysis 

-  Refers  to  a  MIST  case  or  imports  an  archived  MIST  case. 

-  In  optimization,  constraints  can  be  placed  on  flutter  damping. 

Perform  a  state  space  gust  response  analysis  for  a  given  set  of  Mach  number, 
boundary  conditions,  method,  and  model  for  a  particular  gust  spectrum. 

-  Refers  to  a  MIST  case  or  imports  one. 

-  In  optimization,  constraints  can  be  placed  on 


58 


STTR  Final  Rpl050896 


1)  RMS  {u},  {«},  or  {ii}  of  structural  points. 

2)  RMS  section  loads  (forces  and  moments)  of  a  group  of  aerodynamic 
components  such  as  a  wing  with  a  tip  missile. 

3)  RMS  control  command  { or  ( 5}. 

-  Triggers  the  computation  of  {ik)  and  for  the  gust  and  load  modes. 

In  the  ZAERO  and  MIST  cases,  the  user  can  archive  a  pre-defined  collection  of  data  that  will 
enable  the  downstream  (or  dependent)  analysis  cases  to  import  previously  computed  data.  This 
allows  the  output  from  the  computationally  intensive  ZAERO  and  MIST  steps  to  be  used 
repeatedly.  The  MIST  step  will  further  support  an  export  option  to  output  the  ASE  model  to 
MATLAB  or  MATRIX  x. 

The  current  ASTROS  MODES  discipline  will  be  enhanced  to  allow  the  normal  modes  to  be 
archived  for  re-use  as  well.  This  will  preserve  the  normal  modes  for  the  generalized  model 
update.  Similar  to  the  current  ASTROS  capability,  the  use  of  more  than  one  set  of  MODES  in  a 
single  boundary  condition  will  not  be  allowed.  Under  the  proposed  effect,  there  will  be  no 
provision  for  combining  symmetric  and  antisymmetric  modal  models  to  create  an  asymmetric 
system. 

All  other  existing  solution  control  commands  will  be  supported  without  change.  The  existing 
flutter  analysis  features  will  continue  to  be  supported,  but  will  not  make  use  of  any  of  the 
enhancements  presented  in  this  proposal  regarding  controls.  The  feature  will,  however,  be 
modified  to  import  ZAERO  and  modal  data  to  take  advantage  of  those  new  features. 

5.4  MAPOL  Modifications 

The  basic  structure  of  the  ASTROS  MAPOL  sequence  will  be  retained  under  the  proposed 
enhancements.  That  structure  is  shown  in  Fig  (5.4.1)  as  follows: 


PREFACE 


UPDATE 

MODELS 


OPTIMIZATION 


ANALYSIS 


SCREEN  h 


SENSITIVITY 


RESIZE 


ON  CONVERGENCE 


Figure  5.4.1  Basic  Structure  of  the  ASTROS  MAPOL  Sequence. 

One  of  the  proposed  new  disciplines  will  occur  in  the  PREFACE  section  (ZAERO)  while  the 
remaining  four  will  be  supported  within  the  analysis  phase  of  both  the  optimization  and  the 
analysis  segments  (Ref  13).  Support  of  the  MIST  case  within  the  optimization  segment  will 
require  some  special  code  in  MAPOL  to  avoid  re-computing  the  minimum  state  approximation 
^ter  the  first  design  iteration.  Note  that,  if  the  ASEFLUT,  MARGIN,  or  GUST  disciplines  use 
imported  RAA  data,  the  MIST  case  need  not  appear  in  the  OPTIMIZE  packet  of  the  solution 
control.  This  represents  the  anticipated  scenario  in  which  the  ZAERO,  MODES,  and  MIST  steps 
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are  performed  in  the  ANALYSIS  segment  of  an  earlier  ASTROS  execution  and  archived  for  use 
in  optimization. 

The  ASEFLUT,  MARGIN,  and  GUST  disciplines  must  be  coded  differently  in  the 
OPTIMIZATION  segment  than  in  the  ANALYSIS  segment  to  accommodate  the  generation  of  a 
new  set  of  generalized  matrix  data  for  K  and  M  using  the  original  normal  modes.  The  MIST 
output  data  are  assumed  to  be  design  invariant  and  will  not  be  updated  either  due  to  changes  in 
the  plant  model  or  the  controller  model.  In  addition,  the  MAPOL  support  of  the  MODES 
discipline  must  be  enhanced  in  the  optimization  phase  to 

1)  Avoid  recomputing  modes  unless  the  following  disciplines  are  present 

a)  flutter 

b)  modal  frequency 

c)  modal  transient 

d)  a  MODES  case  with  applied  frequency  constraints 

2)  If  the  MODES  discipline  must  be  recomputed,  the  \Khi^  and  {Mhi^  matrices 
should  be  used  rather  than  the  updates  [AT^],  and  [Af^].  TTiis  makes  the  MODES 
discipline  very  inexpensive  on  subsequent  design  iterations. 

The  response  analyses,  ASEFLUT,  MARGIN,  and  GUST  will  be  decomposed  into  a  series  of 
discrete  modules  that  be  called  within  either  the  OPTIMIZATION  or  ANALYSIS  segments. 

Since  the  ASEFLUT  and  GUST  disciplines  share  input,  the  modules  creating  those  matrices  will 
be  separately  addressable  from  MAPOL  to  maximize  code  re-use.  As  part  of  the  proposed  effort, 
a  more  detailed  study  of  the  computational  steps  will  be  performed  to  determine  the  appropriate 
decomposition  for  the  MAPOL  modules  supporting  these  disciplines. 

The  ZAERO  cases,  whether  in  the  OPTIMIZE  or  ANALYZE  packets,  will  be  handled  in  the 
PREFACE.  This  is  similar  to  the  current  ASTROS'  handling  of  DLM  and  CPM  matrix 
generation,  except  that,  under  the  proposed  effort,  the  user  will  have  to  caU  a  ZAERO  discipline 
to  compute  the  unsteady  AIC's  as  opposed  to  the  current  approach  where  the  dependent  analysis 
(e.g.  FLUTTER)  triggers  the  computation  of  the  required  data.  The  current  ASTROS 
methodology  is  incompatible  with  the  ARCHIVE/IMPORT  concepts. 

For  sensitivity  analysis,  the  proposed  enhancements  require  a  substantial  number  of  changes. 
Most  importantly,  the  set  of  design  parameters  has  been  extended  to  include  design  variables  that 
have  no  stiffness  or  mass  sensitivity.  ASTROS  must  be  modified  accordingly  to  only  compute 
the  appropriate  elemental  sensitivities. 

The  new  design  constraint  functions  generated  by  the  ASEFLUT,  MARGIN  and  GUST 
disciplines  must  have  the  appropriate  modules  written  to  compute  the  first  derivatives  of  the 
response  in  the  sensitivity  phase  of  the  optimization  segment  As  for  the  case  of  the  analyses,  a 
study  will  be  performed  to  determine  the  best  module  decomposition  to  use  in  supporting  the 
sensitivities  of  the  related  system  matrices.  Further  changes  in  the  sensitivity  analysis  are 
required  to  support  the  concept  that  the  original  normal  modes  (rather  than  the  current  modes) 
are  the  generalized  coordinates. 

The  DESIGN  module  will  need  to  be  enhanced  to  support  the  new  design  constraints  and  design 
variables.  Also,  the  constraint  screening  modules  ACTCON  and  ABOUND  will  need  to  be 
extensively  modified  to  support  both,  the  new  constraints  and  the  new  design  variables. 
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SECTION  6 


FUTURE  RESEARCH  AND  DEVELOPMENT 

The  success  of  the  Phase  I  contractual  performance  leads  to  the  follow-up  action  of  planning  the 
Phase  n  activity.  Under  the  Air  Force's  STTR  invitation,  a  Phase  n  proposal  is  in  preparation  to 
be  submitted  to  the  Air  Force  shortly.  In  the  Phase  n  proposal,  the  ZONA  Team  proposes  to 
complete  the  research  begun  in  Phase  I  and  fully  develop  the  ZAERO  module  and  the  ASE 
module  in  ASTROS.  Seamless  integration  of  these  modules  into  ASTROS  would  be  the  next 
major  task.  ZONA  believes  that  its  ZAERO  module  would  adequately  support  the 
ASTO.OS/ASE  module  for  a  wide  range  of  applications  in  the  ASTROS  design  and  analysis 
environment.  In  the  course  of  current  and  future  developments  of  ASTROS,  ZONA  intends  to 
play  a  major  role  in  the  enhancement  support  and  maintenance  of  ASTROS/ZAERO  module  and 
cons^uently  the  ASE  module.  Once  the  seamless  integration  of  ASTROS/ZAERO  is  completed 
in  this  phase,  its  commercialization  can  be  carried  out  readily  by  a  tentative  ZONA/UAI 
partnership.  Since  both  parties  have  established  a  customer  base  for  codes  in  ZAERO  and  for 
ASTROS,  a  fruitful  outcome  of  the  initial  commercialization  effort  is  anticipated. 

6.1  Anticipated  Successful  Results  of  Phase  H 

(a)  ASTROS/ZAERO 

i)  Full  development  of  the  ZAERO  module  together  with  the  AGM  module. 

ii)  Seamless-Integration  of  ASTROS/ZAERO  (with  AGM)  will  be  accomplished. 

iii)  Full  documentation  of  ASTROS/ZAERO;  ZONA  will  deliver  four  manuals,  namely, 
the  User's  Manual,  the  Theoretical  Manual,  the  Applications  Manual  and  the 
Programmer's  Manual. 

iv)  Commercialization  of  the  ASTROS/ZAERO  will  be  initiated  by  a  tentative  partnership 
between  ZONA  and  UAI,  following  the  incremental  development  of  the  ZAERO 
model  proposed  in  Phase  n. 

(b)  ASTROS/ASE 

i)  Completion  of  the  research  phase  of  the  ASTROS/ASE  development.  Data  retrieval 
(UAIC/USDA)  by  ASE  module  from  ZAERO  module. 

ii)  ASTROS/ASE  pilot  software  created  and  tested.  Proposed  scenario  cases  and 
industry-supplied  realistic  cases  will  be  tested  and  demonstrated  for  the  capability  of 
ASTROS/ASE. 

iii)  Full  documentation  of  ASTROS/ASE.  ZONA  will  deliver  four  manuals  (same  as 
those  in  section  (a)  part  (iii)  above). 

Under  the  Phase  II  support,  ZONA  realizes  that  the  provided  funding  could  only  be  budgeted  to 
accomplish  the  preparation  for  commercialization  of  ASTROS/ZAERO.  There  appears  to  be 
inadequate  budget  for  performing  the  Seamless-Integration  of  ASTROS/ASE,  thus  prohibiting  its 
immediate  commercialization  after  Phase  n.  However,  ZONA  believes  that  funding  from  the 
private  sector  can  be  obtained  from  the  aerospace  industry  under  Phase  HI  solicitation.  The 
completion  of  the  ASE  as  a  documented  code  and  its  commercialization  will  be  discussed  next. 

6.2  Plans  for  Phase  HI 

(a)  Full  commercialization  of  ASTROS/ZAERO:  Further  improvement  of  the  commercialized 
packaging  of  ASTROS/.ZAERO,  which  includes  the  consolidation  of  software  licensing, 
maintenance  user  serve  policies  and  the  refinement  of  the  documentation  package. 

(b)  Enhance  co-marketing  activities  between  ZONA  and  UAI. 
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(c)  Seek  private  funding  support  from  aerospace  industry  for  further  development  of 
ASTROS/ASE: 

i)  Seamless  Integration  of  ASTROS/ASE. 

ii)  Link  up  with  the  aerospace  industry  for  ASE  applications  to  ongoing  projects  (e.g. 
AAW,  HSCT,  JAST,  etc.) 

iii)  In  this  regard,  ZONA  has  obtained  a  number  of  supportive  letters  on  the 
ASTROS/ASE  development  from  potential  users.  The  likelihood  of  funding  support  is 
evidenced  by  these  responses  from  the  aerospace  community. 

6.3  Significance  and  Impact  of  Phase  II  on  Phase  m 

a)  The  successful  results  of  Phase  n  will  enhance  the  confidence  from  potential  users  in  the 
aerospace  industry.  For  example,  Transonic/Hypersonic  flutter  results  of  Z^RO  and  the 
scenario  and  realistic  case  demonstration  of  the  ASE  capability  could  definitely  serve  this 
purpose. 

b)  The  research  and  development  effort  in  Phase  II  is  aimed  at  setting  a  firm  foundation 
towards  a  unified  methodology  for  aeroservoelasticity  design/analysis,  thus  leading  to  a 
industrial-standard  production  software  which  would  serve  the  aerospace  community 
worldwide. 
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SECTION  7 


POTENTIAL  APPLICATIONS 

7.1  Potential  Technical  Applications 

The  potential  areas  of  applications  are  many  for  ASTROS/ASE.  It  will  provide  the  aerospace 

industry  an  excellent  design/analysis  tool  for  aircraft  such  as  the  High  Speed  Civil  Transport 

(HSCT)  or  the  Active  Aeroelastic  Wing  (AAW),  among  other  ongoing  and  future  projects. 

Other  gross  potentials  of  ASTROS/ASE  include  its  integration  with  Probabilistic  Design 

Methods,  Designing  Smart  Structures  and  its  further  inclusion  of  a  emerging  robust  control 

technology. 

7.2  Potential  Commercial  Application 

a)  As  discussed  in  an  earlier  section  (D.2),  ZONA  Technology  has  enjoyed  the  commercial 
success  of  ZONA51  since  its  release  in  1985.  Now  with  over  35  users  worldwide,  the 
ZONA51  code  has  been  marketed  independently  by  ZONA  since  1985  and  jointly 
marketed  by  ZONA  and  the  MacNeal  Schwendler  Corporation  (MSC)  under 
MSC/NAST]^N  (as  Aero  Option  H)  since  1991. 

b)  The  ZAERO  Module 

As  described  previously  (see  Fig  2.1.1),  the  ZAERO  module  contains  four  major  unsteady 
aerodynamic  codes  covering  the  complete  flight  regime  throughout  the  full  Mach  number 
range,  a  capability  surpassing  that  of  the  current  Aero  modules  in  MSC/NASTRAN  and  in 
ASTROS.  Specifically,  ZTAIC  and  ZONA51U  are  the  transonic  and  unified 
hypersonic/supersonic  generalization  of  the  lifting  surface  code  ZONA51,  respectively. 
ZONA6  and  ZONA7  generalize  the  codes  by  inclu^g  geometric  complexity  to  wing-body 
combinations  that  can  represent  realistic  aircraft  configurations  for  aeroelastic  applications. 
Further,  the  UAIC  feature  of  the  ZAERO  allows  it  to  be  readily  interfaced  with  the  FEM 
model  created  by  ASTROS.  A  market  survey  quickly  reveals  that  none  of  the  above 
features  in  terms  of  flow  and  geometric  generalization  or  a  unified  AIC  is  available  in  any 
existing  software.  On  the  other  hand,  ZONA's  market  survey  from  responses  among  our 
ZONA  code  users  and  others  in  the  Aerospace  community  indicates  that  a  software 
package  such  as  the  ZAERO  module  is  lacking  presently  and  would  be  in  high  demand. 

The  built-in  AGM  module  with  ZAERO  will  provide  any  selected  panel  method  a 
convenient  entry  to  ASTROS.  In  this  way,  the  flexibility  of  ASTROS  as  opposed  to  that  of 
NASTRAN  would  be  enhanced  greatly,  thereby  generating  more  users  for  ASTROS. 

c)  ASTROS/ZAERO/ASE 

The  Phase  n  products,  ASTROS/ZAERO  (seamlessly  integrated)  and  ASTROS/ASE  (not 
seamlessly  integrated),  should  stand  out  as  two  unique  software  products  that  are  cutting- 
edge  technology  for  aeroelastic  design/analysis  and  ASE  design/analysis. 

d)  The  commercialization  of  ASTROS/ZAERO  should  resort  to  the  code-licensing  means;  it 
could  be  an  annual  license  or  a  one-time  paid-up  license.  ZONA  has  had  over  ten  years  of 
experience  in  code  licensing  and  UAI  even  more.  The  CREDA  agreement  between 
UAI/AF  should  expand  further  the  worldwide  licensing  of  the  ASTROS/ZONA. 
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7.3  Potential  Use  bv  the  Federal  Government 


A  fraction  of  the  ZONA51  and  ASTROS  users  belong  to  organizations  in  DoD  and  NASA 
(ZONA51  is  used  frequently  by  WPAFB  and  NASA-Langley).  Aeroelasticians  at  Wright 
Laboratory  have  repeatedly  shown  strong  interest  in  the  ZONA  codes  (such  as  ZONA51  or 
ZONAIR)  to  assist  them  in  their  projects  (e.g.  the  AAW  project). 

Since  the  AF  is  the  sponsor  of  Phase  H,  ZONA  intends  to  make  available  the  executables  of  the 
ZAERO  module  to  the  appropriate  divisions  of  WPAFB  for  their  exclusive  usage  on  a  royalty 
free  basis.  No  further  rights  of  the  ZAERO  module  will  be  transferred  to  the  DoD/Federal 
Government. 
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APPENDIX  A 


Contractual  Performance  Report 
of  STTR  AF95T009: 

(Presented  to  Wright  Laboratory  Personnel 
on  March  11-12, 1996,  WPAFB,  Ohio) 
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PROGRESS  OF  PHASE  I 


•  ZONA51U 

-  NASP  Demonstrator  Wing 

•  ZTAIC 

-  Lann  Wing 
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-  Modeled  F-16  Wing 
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TtCHNOLOGY 


AERODYNAMIC  MODEL  OF 
NASP  DEMONSTRATOR  WING 
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V-g  Plot  of  NASP  Demonstrator  at  M  =  5.0 

8  Modes,  Vf  =  2342  KEAS,  Q^=  129  psi,  13  Hz  on  Mode  4 
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V-g  Plot  of  NASP  Demonstrator  at  M  =  10.0 

8  Modes,  Vf  =  2726  KEAS,  =  175  psi,  cOj  =  1 1 .6  Hz  on  Mode  4 
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Minimum-State  (MIST)  Fit  of  GAF  Root-Locus  Plots  of  NASP  Demonstrator 

NASP  Demonstrator,  M  =  15.0  m  -  s  n 


Plots  taken  from:  Scott,  R,C.  and  Potouky,  A.S.,  “A  Method  of  Predicting  Quasi-Steady  Aerodynamics  for 
Flutter  Analysis  of  High  SjKcd  Vehicles  Using  Steady  CFD  Calculations,"  AIAA-93-1364-CP,  1993. 


Root-Locus  Plot  of  NASP  Demonstrator  at  M  =  5.0 


8  Modes,  Q^=  130  psi,  co^=  82  r/s  on  Mode  4 
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Root-Locus  Plot  of  NASP  Demonstrator  at  M  =  10.0 


8  Modes,  Q^=  184  psi,  00^=  71  r/s  on  Mode  4 


REAL(S) 


feZOWA  TECHNOLOeY 

ZTiaiSTTRiKOMIQiS95 

73 


Comparison  of  Flutter  Results 


NASP  Wing 


Aerodynamic 

Method 

Mach  Number 

5  10  15 

qp  ©f  h 

psi  r/s  (Kft) 

qf  cOf  h 

psi  r/s  (Kft) 

qf  h 

psi  r/s  (Kft) 

Piston  Theory 

129  78  18 

184  78  42 

250  72  51 

QSCFD  2d 

169  80  11 

331  81  28 

982  224  22 

QSCFD  3d 

. 

330  82  28 

981  224  22 

ZONA51U 

-  V-g  Method 

-  S-Domain 

129  82  18 

130  82  18 

175  73  42 
184  71  41 

206  70  55 
213  68  54 

h  —  Approximate  matchpoint  altitude. 


Piston  Theory,  QSCFD  2d  and  QSCFD  3d  data  taken  from: 

Scott,  R.C.  and  Pototzky,  A.S.,  "A  Method  of  Predicting  Quasi-Steady  Aerodynamics  for  Flutter  Analysis 
of  High  Speed  Vehicles  Using  Steady  CFD  Calculations,"  AIAA-93-i364-CP,  1993. 
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Root-Locus  Plot  of  NASP  Demonstrator  at  M  =  15.0 


8  Modes,  Q^=  213  psi,  co  68  r/s  on  Mode  4 
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THE  LESSING  WING 


NASA  TND-344  /  H.C.  Lessing  &  J,L.  Troutman  /  M=0.24~1.3 
AR  =  3.0, 5%  Parabolic  Arc 


N 


“T 

.90 


LOCOtion  of  OrifiCn 

%  Chord 

Rod 

orlforj 

5 

75 

10 

12  5 

.20 

22  5 

30 

32.5 

40 

42.5 

50 

52.5 

60 

625 

70 

723 

80 

823 

Rod  orlfces  cre.50  from 
su/'loce  d 


AJJ  d^nensjons  « iocl^ 


ZrU/FQFC-Oct/9$ 


Section 
n  =  0.5 


0.7  0.9 


klOiA  nCHNDLDGY 


The  Lessing  Wing  at  M=0.9 

Phase  Angle  (in  degrees)  of  First  Bending  at  k  =  0.13 


Section  1  Section  2  Section  3 

Ti  =  0.5  Ti  =  0.7  Ti  =  0.9 
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LANN  WING*  CONFIGURATION 


ArWAL-TR-83-3050/A.  STEIGINGA  &  R.  HOUWINK/May  1983 
AFWAL-TR-83-3006/J.D,  MALONE  &  S.Y.  RUO/Fcb  1983 


Lann  Wing  at  M=0.82  Lann  Wing  at  M=0.82 


Lann  Wing  at  M=0.82 _ 

Wing  Pitching  at  62%  Root  Chord,  k=0.205  &  AOA=0.6'’ 


AGARD  STANDARD  445.6  WING* 


TABIZ  I 

MEASnm)  MDAL  TJtBQUiiJmnJ  ASD  PAHEI.  KAS3 


Model  deaeriptloa 

PretpieQCT'y  epa 

Panel  auua, 
alngs 

Panel  apen, 
rt 

Mounting 

Structure 

Model 

fh,l 

^h,2 

^t,l 

^t,2 

'a 

a 

1.250 

v&n 

Solid 

1 

30.lt0 

146.00 

99-10 

243.00 

90.77 

0.02290 

2.500 

Vail 

Solid 

1 

l»k.60 

67-30 

47.70 

117.00 

47.54 

.14347 

2 

14.10 

69.30 

50.70 

127.10 

50.68 

-14650 

2,500 

V&U 

Vealened 

1 

9-70 

47.00 

55.00 

89.50 

34.69 

.13758 

2 

10.10 

49.00 

35-00 

89.00 

35.69 

.13565 

3 

9-60 

50.70 

38.  ID 

90.50 

38.09 

.12764 

>1. 

9-70 

51.00 

36.00 

96.50 

37.88 

.12764 

5 

9.60 

54.20 

39.20 

96.50 

39.07 

.12143 

6 

10.00 

51.20 

30.90 

98.50 

30.77 

.13826 

3.750 

Vail 

Solid 

1 

6.60 

40.90 

32-30 

79-70 

32.29 

.51304 

1.167 

Sting 

Solid 

Left 

30.60 

143.00 

99.00 

268.00 

96-67 

•.01899 

Right 

51.20 

143.00 

97.00 

251.00 

96.97 

•.01899 

^Calculated  panel  maas. 


Modal  daacrlptloa 

tt 

2ba, 

ft 

2bt, 

ft 

A«a, 

Stnetare 

Mo  eat iBg 

Solid 

Vail 

1.250 

0.917 

a.60a 

*3.15 

Solid 

Vail 

2.500 

1.833 

1.208 

*3.15 

Vaakaaad 

Vail 

2.500 

1,831 

1.208 

*3.15 

Solid 

Wall 

1.750 

2.750  1 

1.812 

*3.15 

Solid 

Sting 

I.U7 

.856 

.567 

*3.15 

Figure  1.-  Vlng  panel  dimensions. 


*  AGARD  Report  No.  765,  AGARD  Standard  Acroclaslic  Confignmiinn  for  Dynamic  Response  I-Wing  445.6 
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ACp  of  445.6  Weakened  Wing 

M=0.95,  k=0.17,  8.22%  Semi-Span  Station 

-O.IS 
•0.1 
•0.05 
C^o  0 
o.os 
0.1 
O.IS 

0  0.2  0.4  0.6  0.8  1 

xJc 

Steady  Cp  Computed  by  CAPTSD 


Mode  1  Mode  2 
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ACp  of  445.6  Weakened  Wing 

M=0.95,  k=0.17,  52.45%  Semi-Span  Station 


0  0.2  0.4  0.6  0.3  1 

x/c 


Steady  Cp  Computed  by  CAPTSD 
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I^ZOIWTECHNOIOEY 


Minimum-State  (MIST)  Fit  of  GAF 
445.6  Weakened  Wing,  M  =  0.9 

ZTAIC  Result 


445.6  Solid  Wing  Flutter  Results 


est 
Cases 


(slug/ft3)  (Hz) 


0.00357 


0.00320 


Wind  Tunnel 
Data* 


COf  Vf 

(Hz)  (ft/scc) 


ZONA6 

(Linear) 


©f  Vf 

(Hz)  (ft/sec) 


CAPTSD* 

(Nonlinear) 


©f 

(Hz) 


25.8  435.0 


26.2  472.1 


Minimum-State  (MIST)  Fit  of  GAF 
445.6  Weakened  Wing,  M  =  0.95 

ZTAIC  Result 


Aieu.6euii  AjbuiSeuji 


Root-Locus  Plot  of  445.6  Weakened  Wing  at  M  =  0.9 


4  Modes 


REAL(S)  REAL(S) 
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IMAG(S)  (rad/s) 


Root-Locus  Plot  of  445.6  Weakened  Wing  at  M=0.95 


ACp  of  Modeled  F-16  Wing 

M=0.925,  k=1.0,  5%  Semi-Span  Station 


«/c 


Steady  Cp  Computed  by  CAPTSD 

Mode  1  Mode  2 
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ACp  of  Modeled  F-16  Wing 

M=0.925,  k=1.0, 55%  Semi-Span  Station 


x/c 


Steady  Cp  Computed  by  CAPTSD 

Model  Mode  2 
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ACp  of  Modeled  F-16  Wing 

M=0.925,  k=1.0. 95%  Semi-Span  Station 


Steady  Cp  Computed  by  CAPTSD 

Mode  1  Mode  2 
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ZONA51U  /  ZTAIC  ARE  FOR  ASTROS 


•  AIC  Formulation 

-  Ideal  for  MDO 

•  k-Domain  Based  Aerodynamics 

-  For  All  Conventional  Flutter  Methods 

•  S-Domain  Aerodynamics  For  ASE 

-  Through  Minimum  State  Technique 
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ZONA  AERODYNAMIC  MODULE:  UAIC  —  VSDA 


ZONA6 

0/ 


ZONA 

AERO 

MODULE 


|ZONA7 


0 


ZONA51U 


— ( — I — ; — 

5 

MACH  NUMBER 


-A 

10 


DLM 

\h 


Msa 

NASTRAN 

AERO 


NONE 

/ 


ZONA5I 


10 


MACH  NUMBER 


ASTROS 

AERO 

MODULE 


NONE 
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UNIFIED  S-DOMAIN  AERODYNAMICS  (USDA)  FOR 
ASTROS  ENGINEERING  APPLICATION  MODULES 


86 
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UNIVERSAL  AERODYNAMIC  GEOMETRY  MODULE  (AGM)  I 


GOAL:  to  provide  a  universal  set  of  input  definition  for  aerodynamic 
modeling  of  all  panel  methods  such  as: 

-  High-Order  Surface  Panel  Methods:  PANAIR,  QUADPAN,  ZONAIR,  etc. 

-  Lifting  Surface  Methods:  DIM,  ZONA51,  etc. 

-  Hybrid  Panel  Methods:  ZONAIR  (Thin  Wing  Option) 

^ZONAUCHfJOLDSY 
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UNIVERSAL  AERODYNAMIC  GEOMETRY  MODULE  (AGM)  II 


Major  Tasks 


•  Define  New  Bulk  Data  Entries:  AGRID,  ATRIA,  AQUAD, 
AQDMEM,  ABAR  for  surface  panels  and  CAER07  for  lifting 
surface  panels. 


•  The  New  Bulk  Data  Entries  will  fully  represent  the  Geometry  Data 
required  by  all  panel  methods  such  as  PANAIR,  QUADPAN, 
ZONAIR,  etc. 
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ZONA  PRODUCTS 


ziuinsoiim 


•  Z0NA51 

•  Z0NA51U 

•  Z0NA7 

•  Z0NA6 

•  ZTAIC 

.  ZONAIR 

•  ZSTORE 

^ZOfiATlCHTJOLOGY 


ZONA  PRODUCTS 


-  ZONA51  code  generates  unsteady  supersonic  aerodynamics  for 
lifting  surfaces. 

-  ZONA51 U  code  generates  unified  supersonic/hypersonic 
aerodynamics  for  lifting  surfaces. 

-  ZONA7IZONA6  codes  generate  unsteady  supersonic/subsonic 
aerodynamics  for  aircraft  configurations  with  external  stores. 

-  ZTAIC  code  generates  unsteady  transonic  aerodynamics  and 
transonic  AIC's  for  wing  planforms  and  lifting  surfaces. 

-  ZONAIR  code  a  high  order  panel  code  for  computing  whole 
aircraft  aerodynamics  with  external  stores  including  structural 
flexibility  effects  at  subsonic  and  supersonic  speeds. 

-  ZSTORE  code  generates  store/aircraft  aerodynamics  and  predicts 
store  trajectory  during  the  ejection  and  separation  phases. 
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Nonlinear 

Transonics 


Various  Limits  for  Hypersonic/Supersonic  Flow 


^ZOfUHCHNOLDGY 


Damping-in-Pitch  Derivative  vs  Thickness 


Diamond  Profile  ( x  =  tan  G  ) 


M  =  10.0 
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Unsteady  Newtonian 


NLR  WIND  TUNNEL  MODFJ, 

WING  +  TIPTANK  +  PYLON  +  STORE 


MKr  ter.  PLAMC 


MCA 

6SA^OOH,$  MOOinCP 


tSja  CHO«P 


ViKf 


LMSAi;  pJMrH5lC>JJ  IH 


NLR  WING  TIPTANK  CONFIGURATION 
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NLR  Wing+Tiptank 

M=0.45,  k=0.305,  Xq=0.15Cr 


Present  =  Z0NA6 


Unsteady  Cp  Along  Tiptank  at  9  =  202.5’ 
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NLR  WING+TIPTANK+PYLON+STORE 
CONFIGURATION 


0  ZOHATECHHOLOCY 


AERODYNAMIC  MODELING 
WING  BOXES  =  90 
PYLON  BOXES  =  24 
TANK  PANELS  =  264 
STORE  PANELS  =  216 
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NLR  Wing+Tiptank+Pylon+Store 


M=0.45,  k=0.305,  x„=0.15Cb 


Present  =  ZONA6 


Unsteady  Normal  Load  (C^ )  Along  Store 


NLR  Wing+Tiptank+Pylon+Store 


Unsteady  Cp  Along  Store,  M=0.45,  k=0.305,  Xo=0.15Cj^ 
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NLR  WING  WITH  TTP-MT.SSn,F, 

Spanwise  Airload  on  Wing,  M=0.€,-K=0.2 


Development  of  ASTROS/ASE  Module 

-  Development  of  Software  Design  Blue-print:  Phase  I 

-  Submodule  Development/Total  Module  Integration:  Phase  II 

•  OVERALL  CAPABILITIES 

•  APPROACH 

•  ASE  DESIGN  SCENARIOS 

•  PROGRAM  FLOW  CHART  OF  ASE  MODULE 

•  GROSS  POTENTIALS 

95  ftztmncHNOLOGY 
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OVERALL  CAPABILITIES 


•  Inclusion  of  Multi-Input  and  Multi-Output  (MIMO)  Control  System 

•  Provide  Closed-Loop  Robust  Stability  Analysis  by  Means  of  Gain 
Margins,  Phase  Margins  and  Singular  Values 

•  Add  Continuous  Gust  Response  Capabilities 

•  Stability  and  Gust  Response  Constraints  in  the  Structural  Design 
Optimization 

•  Inclusion  of  User-Defined  Control  Parameters  of  a  Given  Control 
Law  in  the  MDO  Process 

•  Export  an  Efficient  State-Space  Representation  of  the  ASE  System 
to  MATLAB/MATRIX  X  for  Subsequent  Analysis  and  Control 
Synthesis 

The  above  new  features  are  applicable  to  open-loop  as  well  as  closed  loop  systems. 

TECHNOLOGY 
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APPROACH 


•  Rational  Approximation  of  Unsteady  Aerodynamic  Forces 
Provided  by  the  ZONA  Aero  Module 

-  Minimum  State  Technique 

-  Roger’s  Approximation 

•  State-Space  Formulation  of  the  ASE  System 

•  Control  System  is  Represented  by: 

-  Polynomial  Transfer  Functions 

-  State-Space  Realization 

•  Gust  Filter  is  Defined  such  that  a  White-Noise  Input  Produces 
an  Approximation  of  Either  Dryden's  or  von  Karmen's  PSD 
of  Atmospheric  Continuous  Gusts 

•  Residualization  of  Structural  States  for  Modal-Size-Reduction 

^ZONA  TECHNOLOGY 
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ASE  DESIGN  SCENARIOS 


Three  Design  Scenarios  are  Presented  as  Follows: 


•  AFW  Roll  Performance  with  Flutter  Constraints 

•  Trade-Off  Studies  of  Gust  Load  Alleviation  by 
Passive  and  Active  Means 

•  Store  Flutter  Suppression  with  Structural 
Uncertainties 


^ZOfW  TECHWOLOfiY 
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ASE  DESIGN  SCENARIOS  (I) 


AFW  Roll  Performance  with  Flutter  Constraints 


•  The  AFW  wind  tunnel  model  has  two  leading  edge  &  two  trailing 
edge  control  surfaces  to  meet  the  roll  performance  requirements. 

•  An  initial  control  system  design  satisfies  these  roll  constraints 
but  introduces  antisymmetric  flutter. 

•  With  the  proposed  ASTROS/ASE  module,  the  control  gains  of 
the  four  control  surfaces  are  added  to  the  list  of  design  varaibles. 

•  The  new  ASE  module  and  the  static  aeroelastic  module  are 
employed  to  optimize  the  control  gains  and  the  structural  design 
variables  for  satisfying  the  flutter  and  roll  performance  constraints 
with  a  minimal  weight  penalty. 

^ZCNA  TtCHNOLOSY 

ZrL»jrTtC/WU’us02l496  97  ^ 


ASE  DESIGN  SCENARIOS  (II) 


Trade-Ojf  Study  of  Gust  Load  Alleviation  by  Passive 
and  Active  Means 


•  A  drone  aircraft  designed  by  ASTROS  Version  11  is  experiencing  excessive 
wing  root  RMS  bending  moment  in  open-loop  continuous  gust  analysis. 

•  With  the  new  ASE  module,  the  baseline  structure  is  used  as  a  starting 
point  for  two  design  studies; 

♦  Passive  Design  -  Add  gust  response  constraints  and  repeat  the  structural 

design  process. 

♦  Active  Design 

-  Extract  the  plant  state-space  model  for  designing  a  load  alleviation 
control  system  outside  of  ASTROS. 

-  Due  to  the  control  system  limitations,  the  control  law  helps  but  does 
not  fully  satisfy  the  gust-response  requirements. 

-  Redesign  the  structure  in  the  presence  of  the  control  system  with 
gain-margin,  phase-margin  and  singular-value  constraints  added  to  the 
list  of  constraints. 

•  Supply  data  for  trade-off  studies  of  structural  changes  vs.  complexity  of  the 

control  system.  ^zonatjchnolosy 

ZTtSt^RtWLPritOZ  H96 


ASE  DESIGN  SCENARIOS  (III) 


Store  Flutter  Suppression  with  Structural  Uncertainties 


•  Uncertainties  of  the  stores/pylons  structural  properties  require 
a  robust  control  system  design  for  store  flutter  suppression. 

•  The  state-space  matrices  of  aircraft  with  external  stores  and  the 
derivatives  of  these  matrices  with  respect  to  the  structural 
properties  of  the  stores/pylons  will  be  generated  by  the  new 
ASE  modules. 

•  MATLAB  p-analysis  tool  box  will  then  be  used  to  design  a 
control  system  which  suppresses  flutter  over  a  range  of 
parameter  uncertainties. 

•  The  designed  control  law  will  be  introduced  to  ASTROS  for 
the  verification  with  various  store  models  at  different  flight 
conditions. 
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•  Integration  with  Probabilistic  Design  Methods 

•  Introduction  of  Smart-Structure  Technology 

•  Inclusion  of  Emerging  Robust  Control  Technique 
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